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ABSTRACT 

'this  grant  supported  our  work  on  semiconductor  alloy  theory. 
Many  useful  results  have  been  obtained,  including:  (1) 
generalization  of  Brooks'  formula  for  alloy-scattering 
limited  electron  mobility  to  including  multiple  bands  and 
indirect  gaps,  (2)  calculation  of  SiGe  alloys  band  structure, 
electron  mobility  and  core-exciton  binding  energy  and 
linewidth,  (3)  comprehensive  calculation  of  bond  energy,  bond 
length  and  mixing  enthalpy  for  all  III-V  and  II-VI  pseudo¬ 
binary  alloys,  (4)  development  of  a  statistical  theory  which 
shows  a  non-random  distribution  of  atoms  in  most  alloys,  (5) 
studying  the  sensitivity  of  defect  levels  to  band  structures 
and  impurity  potentials,  (6)  a  study  of  the  dipolar 
contribution  to  the  mixing  energy  and  its  implication  to  the 
long-range  order  in  alloys,  e.g.,  GaAlAs,  and  (7)  a  model 
which  allows  a  simple  but  detailed  calculation  of  alloy 
band-edge  properties. 


ACCrirx.n  O.ia  Ex 


SUMMARY 


This  grant  supported  our  work  on  semiconductor  alloy 
theory.  Through  interactions  with  Dr.  Arden  Sher's  group  at 
SRI  International,  many  useful  results  have  been  obtained  under 
this  grant.  Below  we  summerize  th  highlights.  The  details 
will  be  discussed  in  the  publications  enclosed  with  this 
report . 


A.  Generalized  Brooks'  Formula  and  the  Electron 
Mobility  in  SiGe  Alloy 

Although  Brooks'  formula  has  been  used  widely  for 
calculating  the  alloy-scattering  limited  electron  mobility,  we 
show  that  this  formula  is  only  valid  for  a  direct-gap 
semiconductor.  There  are  also  questions  about  the  scattering 
parameter  and  the  effective  mass.  We  generalized  the  formula 
for  indirect-gap  alloys  with  multiple  bands  and  applied  it  to 
SiGe  alloy.  Our  results,  correlated  well  with  experiments, 
showed  that  the  electron  mobility  drops  fast  with  alloying. 

The  mobility  has  a  dip  at  15%  Si  concentration,  corresponding  a 
transition  from  the  X  to  the  L  edge. 

B.  Bond  Lengths,  Lattice  Relaxation  and  Mixing 
Enthalpies  in  Semiconductor  Alloys 


We  treated  the  problem  with  a  model  which  combines 
Harrison's  bonding  theory  with  a  valence  force  model  and  an 


elastic  continuum.  While  the  local  strain  is  the  main  driving 
force  for  the  bimodal  bond-length  distribution  in  pseudo-binary 
alloys  found  in  the  EXAFS  experiment,  we  found  that  the 
chemical  shifts  arising  from  different  bond  lengths  and 
polarities  of  the  constituent  bonds  can  have  a  significant 
contribution  to  the  mixing  enthalpies.  In  fact,  the  sizable 
negative  values  of  the  chemical  shifts  in  the  cation 
substitutional  alloys,  e.g.  Ga (x) In (1-x) As,  may  be  important 
for  stabilizing  the  mixture.  We  also  deduced  a  simple 
criterion  for  separating  miscible  from  immiscible  alloys. 

C.  Sensitivity  of  Defect  Levels  to  Host  Band 
Structures  and  Impurity  Potentials 

Our  calculation  of  defect  levels  for  more  than  30 
impurities  in  CdTe  showed  that  the  discrepancies  caused  by 
different  host  band  structures  and  impurity  potentials  ranged 
from  less  than  0,1  eV  to  the  whole  band  gap  (1.6  eV) .  This 
result  casts  some  doubt  about  the  quantitative  nature  of  the 
empirical  tight  binding  method  for  deep-level  studies.  We  also 
suggested  ways  to  improve  the  theory. 

D.  SiGe  Alloys  -  Band  Structure  and  Core-Exciton 

SiGe  binary  alloy  has  regained  research  interests  recently, 
because  it  has  a  potential  for  high-speed  devices  in  the 
strained  superlattice  configuration.  We  have  applied  our 
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technique  to  obtain  high-quality  band  structures  for  Si  and  Ge 
and  have  performed  alloy  calculation  for  S i  ( x ) Ge  ( 1-x ) .  The 
results  have  been  checked  against  available  optical  data  and 
have  been  applied  to  the  mobility  calculation  mentioned 
earlier.  Another  interesting  result  is  that  the  alloy  band 
parameters  allowed  us  to  correlate  the  Si  2p  core-exciton 
binding  energy  with  its  linewidth  in  the  alloy.  The  observed 
minimum  in  the  linewidth  near  x=0.15  can  be  explained  as  a 
result  of  competition  between  an  intrinsic  broadening  due  to 
screening  and  an  extrinsic  alloy  scattering.  The  most 
reasonable  binding  energy  in  pure  Si  was  found  to  be  0.15  eV. 
This  work  thus  has  helped  resolve  the  controversy  about  the 
unusually  large  binding  energy  for  the  core  exciton  in  Si. 

E.  Calculation  of  the  Alloy  Band-Edge  Properties 

There  are  two  basic  techniques  for  calculating  the  band 
structure  in  semiconductor  alloys:  empirical  formula  and 
detailed  theory  such  as  the  coherent-potential  approximation 
(CPA).  The  formal  is  handy  to  use,  but  provides  very  little 
insight.  The  latter  is  usually  very  time  consuming,  and  for 
energies  near  band  edges  the  numerical  accuracy  is  often 
questionable.  For  1 1  I -V  and  most  of  the  II-VI  pseudo-binary 
alloys,  these  difficulties  can  be  circumvented  by  using 
perturbation  theory.  Moreover,  the  CPA  results  are  sensitive 


the  band  models  used  for  the  constituent  compunds.  By 
examining  the  overall  structures  of  the  Green's  function, 
partial  densities  of  states,  and  the  tight-binding  interaction 
parameters  involved  in  the  band-edge  states  and  in  the 
perturbation  theory,  we  have  concluded  a  procedure  for  a  simple 
but  detailed  calculation  of  the  various  contributions  to  the 
band  gap  bowing  and  the  low-field  electron  mobility  in  these 
a  1 loys . 

F.  Dipolar  Contribution  to  Alloy  Mixing  Energy 
and  Its  Implication  to  Long  Range  Order 

The  alloy  formation  energy  (F.E.),  defined  as  the 
difference  between  tha  alloy  total  energy  and  the  average  of 
the  constituents'  values,  and  mixing  entropy  are  two  competing 
factors  which  govern  the  alloy  thermodynamics.  While  a 
negative  E  favors  a  long-range  order  (LRO)  at  low  temperature, 
a  positive  F.E.  tends  to  cause  spinodal  decomposition.  It  was 
suggested  recently  that  the  long-range  electrostatic 
interaction  among  cations  and  anions  contributes  to  a  negative 
F.E.,  and  hence  serves  as  a  driving  force  for  the  LRO  found  in 
GaAlAs  grown  by  MBE  and  MOCVD.  We  showed  analytically  that, 
when  all  terms  are  included,  this  dipolar  contribution  to  F.E. 
is  positive.  Therefore,  it  is  not  a  force  to  stabilize  the  LRO 
phase.  However,  the  value  of  the  dipolar  F.E.  in  the  ordered 
phase  is  smaller  than  that  in  the  disordered  one.  Since  the 
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Semiconductor  pseudobinary  alloys:  Bond-length  relaxation  and  mixing  enthalpies 


A.-B.  Chen 
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(Received  1  March  1985) 

Harrison's  bonding  theory,  the  valence  force  field  (VFF),  and  an  elastic  continuum  are  combined 
in  a  study  of  the  substitution  energies  A,  and  local  (first-shell)  bond  lengths  d,  of  isoelectronic  im¬ 
purities  in  semiconductors.  Explicit  expressions  for  A,  and  d|  are  denved,  which  enable  us  to  ab¬ 
sorb  measured  elastic  constants  into  the  calculation  and  to  study  the  chemical  effects  arising  from 
differences  in  the  covalent  radii  and  polarities.  Several  models  based  on  VFF  alone  are  also  derived 
for  comparison.  The  full  theory  and  at  least  five  VFF  models  are  found  to  produce  impurity  bond 
lengths  in  excellent  agreement  with  experiment.  The  substitution  energies  are  shown  to  provide 
good  estimates  of  the  mixing  enthalpies  fl  of  pseudobinary  alloys  and  to  predict  miscibility  gaps 
properly.  The  chemical  shifts  in  Cl  are  found  to  be  negative  for  most  cation  alloys  but  positive  for 
anion  substitutions. 


I.  INTRODUCTION 


TTie  discovery  of  a  bimodal  distribution  of  the  nearest- 
neighbor  bond  lengths'  in  Ga;,In|  .j,As  has  sparked  con¬ 
siderable  interest  in  the  bonding  nature  of  semiconductor 
alloys.^”*  This  finding  has  changed  the  conventional  pic¬ 
ture  of  the  alloy  crystal  bond  configuration,  which  has 
far-reaching  implications  about  the  electronic  structure, 
structural  stability,  and  thermodynamics  of  these  materi¬ 
als.  Because  of  the  complexity  of  both  the  structural  and 
the  potential  disorder  in  these  alloys,  ab  initio  band- 
structure  techniques  have  not  yet  evolved  to  a  stage  suit¬ 
able  for  direct  calculations.  Therefore,  we  have  extended 
Harrison’s  bonding  theory  to  study  the  alloy  structural 
properties.'  *’  In  this  paper,  we  apply  an  intermediate  ver¬ 
sion  of  the  theory  to  the  dilute-limit  case  of  an  isoelect¬ 
ronic  impurity. 

A  particularly  useful  application  of  the  theory  is  its 
pierturbation-expansion  form,  in  which  measured  elastic 
constants  are  incorporated  to  obtain  accurate  results. 
This  form  is  also  useful  for  comparison  with  other  previ¬ 
ously  published  models'  that  are  based  on  the  valence- 
force-field  (VFF)  (Ref.  10)  model  alone.  Thus,  all  the 
factors  influencing  bond-length  relaxation,  e.g.,  strains, 
boundary  conditions,  and  chemical  effects,  can  be  studied. 
The  ability  to  incorporate  the  chemical  effects  is  one  ma¬ 
jor  difference  between  this  theory  and  other  VFF  models. 

The  remainder  of  the  paper  contains  the  fi'llowing  sec¬ 
tions  Sec.  II  descnbcs  a  theory  for  calculating  impurity 
substitution  energies.  Section  III  casts  the  theory  into 
perturbation  form  and  combines  it  with  a  valence  force 
field  and  an  elastic  continuum.  Several  VFF  models  are 
derived  m  Sec.  IV,  The  mtxlifications  of  numerical  re¬ 
sults  due  to  chemical  effects  on  kval  bond  lengths  and  al¬ 
loy  mixing  enthalpies  are  summari/ed  and  discussed  in 
Sec.  V.  Conclusions  are  drawn  in  the  last  section.  Sec.  VT. 
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II.  IMPURITY-SUBSTITUTION  ENERGY 

Consider  the  problem  of  substituting  an  isoelectronic 
atom  A  for  &  B  atom  in  a  zinc-blende  compound  BC  (e.g., 
In  substitutes  for  Ga  in  GaAs,  as  shown  in  Fig.  I).  In 


0  on  0 

V  3  V  3 


no  1  A  sketch  the  flattened  picture  (^f  a  single  impuniy 
.4  in  a  HC  Cvimpound.  The  positions  and  displacements  for 
those  atoms  labeled  are  used  in  .Appendix  H 
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general,  the  bond  lengths  d),  d;,  dj, ...  for  the  first-, 
second-,  and  third-shell  bonds  surrounding  the  impurity 
are  different  from  the  equilibrium  values  of  either  the 
pure  BC  compound,  denoted  as  d,  or  the  "impurity”  com¬ 
pound  AC,  denoted  as  d;  =d(  1  — 8o).  If  A  starts  being  a 
free  atom  and  B  also  ends  being  a  free  atom,  then  the  en¬ 
ergy  difference  between  the  final  and  initial  states  is  de¬ 
fined  as  the  substitution  energy  and  is  given  by 


A. 


(1) 


l^def +  ^8  )  — (£pu„  -t-E^  ) 

^  and  Eg  are  free-atom  energies  for  A  and  B 


where  _ 

respectively,  E^cf  is  the  total  energy  of  the  semiconductor 
with  a  defect  as  sketched  in  Fig.  1,  and  is  that  of  the 
pure  host  BC  crystal.  Equation  (1)  can  be  written  as 

“Eal  +  (£^dis~£purc) 

=  A8— (E^— Egl-t-  Ajji,,  ,  (2) 

where  we  have  added  and  subtracted  a  term  which  is 
the  total  energy  of  a  distorted  BC  compound  with  all  the 
atoms  held  at  the  positions  specified  in  Fig.  1,  except  that 
the  central  atom  is  a  fi  atom.  Clearly,  Ay,  =E(,<.f  —  Ej^  is 
a  replacement  energy,  and  the  distortion  energy 
^dis=  £^<iis  “  ^purc  is  the  energy  required  to  deform  a  pure 
BC  crystal  from  its  equilibnum  lattice  to  that  specified  in 
Fig.  1.  A/,  contains  all  the  chemical  terms  that  arise  from 
different  bond  lengths  and  polarities  between  AC  and  BC 
bonds. 

Ay,  can  be  treated  most  easily  by  Harrison's  bonding 
theory.’  In  this  theory,  the  energy  per  bond  relative  to  the 
vacuum  state  is 


Ej,  =  ICh  +  F()  +  be  *  -1  ()£„ 


(31 


where  E^  is  the  energy  of  the  bonding  state  constructed 
from  the  two  hybrid  orbitals  facing  each  other  along  the 
bond  direction 


E»  =  jlE^-t-E^'l-IFt  +  ki 


2,1/2 


14) 


with  Ey,'"  and  eJ  being  the  energies  of  the  anion  and  cation 
hybnd  orbitals,  respectively.  The  antibonding  energy  e^ 
has  the  same  form  as  in  Eq.  i4),  except  with  a  plus  sign. 
I'l  is  culled  the  covalent  energy,  which  is  the  total  elec¬ 
tronic  Hamiltonian  matrix  element  between  the  two  hy¬ 
brids  in  question,  and  the  polar  energy  F,  is  the  differ¬ 
ence  F-,  =  7 1  eJ  -  Ey,^ '.  The  E*  and  c,;,,  the  metallization 
energies,  are  the  shifts  of  the  b<inding  level  caused  by  in¬ 
teractions  with  the  neighboring  antibonding  states,  where 
-I-  and  -  indicate  whether  the  common  adjacent  atom  is 
a  cation  or  anion.  For  example,  this  term  for  an  AC  bond 
labeled  by  d,  in  Fig.  1  due  to  an  antibonding  state  labeled 
by  d-  is  given  by 

1  •  ' 


,-l2.1  )  =  - 


£6I1'-e;(2) 


(5) 


where  .4  and  B  denote  AC  and  BC  bonds,  respcctivelv. 

=  x'eJ  -£p'.  with  c‘  and  i^  being  the  s-  and  p-term 
values  of  the  common  adjacent  atom  C.  (V(l)  is  the 
probability  amplitude  of  finding  an  electron  in  the  hybnd 
orbital  of  the  C  atom  in  the  bonding  state  of  an  /4Cbond 


with  a  bond  length  d,.  whereas  t/*(2)  is  the  correspond¬ 
ing  probability  amplitude  for  the  antibonding  state  of  a 
BC  bond  of  bond  length  di-  Finally,  Vq  is  a  repulsive 
pair  potential  required  to  prevent  the  crystal  from  collaps¬ 
ing  and  to  guarantee  a  correct  equilibrium  bond  length. 
The  local  perturbation,  Eq.  (5),  is  applicable  because  the 
square  root  of  the  numerator  is  much  smaller  than  the 
separation  between  the  antibonding  and  bonding  levels 
and  the  valence  band  is  completely  filled,  so  the  interac¬ 
tion  between  the  bonding  states  only  spreads  the  Ey,  levels 
into  bands  without  affecting  the  center  of  gravity  of  the 
occupied  states. 

The  replacement  energy  Ay,  of  Eq.  (2)  can  now  be  writ¬ 
ten  explicitly: 

Ay,  =  4[  2e^  1 )  +  FiJ  ( 1 1  +  6E;i;l  (1.1)  +  6e»vJ.(  2, 1 ) 

-  2Ef ( 1 )  -  F^(  1 )  -  6e®' *8 ( 1 . 1 )  -  6e«' ( 2, 1 ) 

+  6£;i/<l,2)-6E*'*r(1.2)]  .  (6) 

The  distortion  energy  Aj^  of  Eq.  (2)  now  involves  only  BC 
bonds  of  different  bond  lengths.  It  can  be  treated  with  ex¬ 
actly  the  same  procedure  for  any  given  set  of  bond-length 
distributions.  Thus,  a  straightforward  energy  minimiza¬ 
tion  procedure  can  be  carried  out.  The  accuracy  of  this 
procedure,  however,  depends  in  turn  on  the  accuracy  of 
scaling  rules  for  Fj  and  Fq  and  the  input  parameters.  At 
present,  Harrison's  model’  with  Fy  a  l/d‘  and  F^a  l/d* 
and  his  universal  parameters  are  only  .semiquantitative. 
We  are  improving  the  quantitative  nature  of  the  theory  so 
that  the  full  theory  will  yield  accurate  predictions  of  the 
structural  and  thermal  properties  of  semiconductor  de¬ 
fects  and  alloys. 


III.  PERTURBATION  EXPANSION. 

valence  force  field, 

AND  ELASTIC  MEDIUM 

As  pointed  out  earlier,  a  perturbation  expansion  of  the 
theory  is  instructive.  This  is  feasible  because  the  differ¬ 
ences  AF,  =  F,(/1C)-F2(BC)  and  AF,=  F,(/1C) 
—  F)(BC)  are  small  compared  to  each  individual  value  for 
many  of  the  isoelectronic  impurities  in  IIl-lV  and  II-VI 


compounds.  To  this 

end. 

Eq.  16)  is  rew  r 

4Ay,  =  E„'(l)- 

Eyf(  1 

)  +  6[e*V<-(2.11 

-t>[E,i:f 

(2.1 1 

-£*+-(1.11] 

(1.21 

-£^..*'2,2.] 

-b[E«| 

1 1,2) 

-£lr(2.2)]  . 

C:c(  1 


.It] 


i7) 

where  E^’i  1 '  and  £yf(  1 )  are  energies  per  bond  in  Eq.  (3) 
for  .4C  and  BC  compounds,  respectively,  with  the  relaxed 
bond  length  dj  =di  1-  6'.  The  difference  between  these 
energies  Ey,  and  the  corresponding  values  at  their  respec¬ 
tive  equilibnum  bond  lengths  d,  and  d  are  just  the  strain 
energies  per  Ixind  in  uniform  deformation; 


E„’(  1 1 


Eyf  I  1 '  - 


.  I  2\  }Bjd,id]  d. 


Eh  d>^ 


}Bdui: 


“I 

.%  Vv 


■  f 

'  ‘A 

• I-. I 


v'  *v 
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f?  f  I 


s'  V- 


9  I 
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■  m 


l  ^  -Ok  , 
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where  B/  and  B  are  the  bulk  moduli  for  the  impurity  AC 
and  host  BC  crystals.  The  rest  of  the  terms  in  Eq.  (7)  are 
all  due  to  changes  in  t„  caused  by  the  differences  AKj 
and  A3'2-  shall  use  Harrison's  scaling  rules  to  deduce 
them.^  Expanding  Eq.  (7)  to  second  order  in  AT}  and 
Ad=d2— d|,  we  write 

6(e*l(2.1)-ei;^(l.l)I 

=/;  Ad  -g/AK3  +  /i,(Ad)2-  W,  Ad  AKj  + 1/,(  AFj)^  , 

(9) 

where  fj,  gj,  and  so  on,  are  appropriate  derivatives 
evaluated  for  the  impurity  crystal  AC.  When  similar  ex¬ 
pansions  are  made  for  the  rest  of  the  terms  in  Eq.  (7),  it 
becomes  [with  di=d(l— 8),  d/=d(l— Sq)] 

T  A*  =  AE*  -I-  (/,  -/)Ad -^g,-g)^V^+{h,+h){^^d)^ 

-^w,  +  w)^d^Vi+^u,  +  u){^yi)^ 

-(-2t/3B/d^(6-8o)^-2t/3Bd^8^  ,  (10) 

where 

A£*=£fc'*(d,)-£:4*(d)-i(e^-efl)  (11) 

is  just  the  difference  in  the  binding  energy  per  bond  be¬ 
tween  the  "BC'  and  “AC’  crystals.  In  Eq.  (10)  the  coef¬ 
ficients  f,g  without  a  subscript  are  those  for  the  host  BC 
system.  It  is  convenient  to  define  an  excess  energy 
A£  =  Aj/4  — Af*.  which  is  the  extra  energy  per  bond  re¬ 
quired  for  the  impurity  substitution  over  and  above  the 
binding-energy  difference  between  the  BC  and  AC  crys¬ 
tals.  The  binding-energy  difference  accounts  for  much  of 
the  substitution  energy;  however,  the  correction  measured 
by  the  excess  energy  can  be  significant.  The  excess  energy 
results  from  strain  energies  and  chemically  driven  charge 
redistributions  around  the  defect.  Using  ^s.  (2)  and  (10) 
and  defining  F=fi—f  and  G=g[—g,  we  can  write  A£ 
up  to  second  order  in  A  and  Ad  as 

A£  =  2t/'3B,d ’(8 -6o)^  -  2t/3Bd^6^ -(- £  Ad  -  G  A  k-j 

-(-//(Ad)^-(-Ik'AdAI'3-(-t/(AK3)^+-J-Ad„  ,  (12) 

v/hert  H  =  hi+h,  H'=wi+w,ind  U  —  u,+u. 

To  treat  the  distortion  energy  A^,,,  we  divide  the  crystal 
into  two  regions.  Inside  a  sphere  of  some  radius  R  mea¬ 
sured  from  the  impurity,  the  strain  energy  is  taken  to  be 
the  valence-force-field'®  value: 

+  AS'  2  AylA'd.-d;)]',  (13) 

Bfl  ,  ! 

ly  <  1 1 

where  i  sums  over  all  the  bonds  inside  R  and  the  pairs  in 
the  B  terms  include  those  that  have  adjacent  atoms  inside 
R  and  on  the  boundary.  The  parameters  a  and  B  are 
force  constants  to  be  considered  later  A(dd^)  =  d,d^ 
— dl®''d|®  measures  the  change  of  the  dot  product  be¬ 
tween  bond  vectors  due  to  distortions.  Outside  R  we  as¬ 
sume  an  elastic  continuum  with  radial  displacements  in¬ 


versely  proportional  to  the  square  of  the  distance  from  the 
center.  The  elastic  energy  in  this  medium  can  be  shown 
to  be  (see  Appendix  A) 

A;S“’=£Cw^  (14) 

where  the  effective  shear  coefficient  is 

C =4)r[0.4(C, ,  -  C,2 ) -)- 1 . 2C4« ) 

and  u  is  the  displacement  at  R.  In  view  of  the  fact  that 
the  bonds  d|  and  di  are  coupled  through  the  chemical 
terms  in  Eq.  (10),  the  smallest  logical  radius  R  is  the 
second-shell  atomic  distance,  namely  £=2v^/t/3. 
Atoms  on  this  boundary  have  displacements  of  the  forms 

u=d(y,y.0)/t/3 .  Thus,  u=vl.yd/'/i  and  the 

elastic  energy  in  the  continuum  is 

AiJJ“’=-f  t/ITSCyV’ .  (15) 

The  distortion  energy  represented  by  Eq.  (14)  contains 
six  different  contributions  (see  Appendix  B):  the  bond¬ 
stretching  energy  of  the  four  first-shell  bonds  6a8^d^,  the 
B  terms  from  the  first-shell  bonds,  B^^d^,  the  a  terms 
from  the  second-shell  bonds,  2a(8-(-2y  )^d^,  the  B  terms 
between  the  first-  and  second-shell  bonds,  2P(8-)-y)^d^ 
the  B  terms  among  the  second-shell  bonds,  IBB^d^,  and  fi¬ 
nally  the  B  terms  between  the  second-shell  bonds  and 
those  in  the  continuum, 

||5d^[(38+X,y)2-f  (8-t-A.2y)^)  , 

where  A.,  =40)/2/(  19v/T9)  and  A.2  =  2-40t/2/(  1 1  vTl ). 

To  assemble  all  the  contributions  to  Eq.  (121,  we  need  to 
consider  the  assignments  of  the  elastic  constants  and  the 
force  constants  a  and  B  in  VFF.  While  the  experimental 
values"  of  C|,>  C12,  and  C44  can  be  used  for  the  elastic 
constants,  a  and  B  have  to  be  deduced.  If  Martin's  origi¬ 
nal  procedure'®  (also  followed  by  Martins  and  Zunger^)  is 
used,  then  Eq.  (13)  alone  will  not  produce  the  correct  (ex¬ 
perimental)  bulk  moduli.  There  are  small  corrections  due 
to  Madelung  terms,  which  are  hard  to  treat  in  the  case  of 
nonuniform  distortions.  A  simpler  procedure  is  adopted 
here.  We  use  the  experimental  bulk  moduli  for  B,  and  B 
in  Eq.  (12)  and  experimental  elastic  constants  to  calculate 
C  of  Eq.  (14)  and  then  force  a  and  B  Ihe  VFF  to  pro¬ 
duce  the  correct  bulk  moduli  B  and  shear  coefficients 
C|i— C’12.  Such  an  approach  is  also  consistent  with 
Hamson's  bonding  theory’  and  other  approaches  in  which 
the  Coulomb  forces  are  automatically  incorporated  in  the 
band  and  bond  energies,  and  do  not  need  to  be  redundant¬ 
ly  treated.  With  our  procedure,  the  bulk  modulus  is  sim¬ 
ply  related  to  the  force  constants  by  B  =  (3a-(-0)/(4v^), 
Table  1  lists  our  a  and  B  values.  We  want  to  point  out  in 
advance,  however,  that  the  numerical  results  deduced 
from  our  sets  and  those  of  Martin  of  a  and  B  do  not  in¬ 
troduce  differences  more  than  the  present  experimental 
uncertainties  in  the  local  bond  length  (-0.01  A)  and 
the  mixing  enthalpies  (  >0.5  kcal/mole). 

Using  the  above  procedure  and  adding  all  contributions, 
the  excess  energy  per  bond  from  Eqs.  (10)  and  (12)  is  the 
full  perturbation  theory  (FPT)  result 
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TABLE  1.  The  btind  lengths  d  (in  A),  valence  force  constants  a  and  0  (N/m),  shear  coefficients  C  of 
the  continuum  (in  10"  ergs/cm'),  melting  temperatures  T„  (K),  and  Liedcrmann  ratios  l'«  for  the  com¬ 
pounds  used  in  this  paper. 


Compound 

d 

a 

0 

C 

T„* 

Xm 

AlP 

2.367 

44.323 

8.068 

122.396 

1773 

0.070 

AlAs 

2.451 

40.849 

8.717 

112.695 

1873 

0.0’ 3 

AlSb 

2.656 

34.073 

6.900 

85.351 

1323 

0.062 

GaP 

2.360 

44.764 

10.737 

145.921 

1510 

0.064 

GaAs 

2.448 

39.235 

9.159 

121.844 

1738 

0.071 

GaSb 

2.640 

31.876 

7.347 

89.372 

985 

0.055 

InP 

2.541 

40.363 

6.543 

91.785 

1343 

0.059 

InAs 

2.622 

33.203 

5.752 

78.816 

1215 

0.061 

InSb 

2.805 

28.557 

4.891 

60.721 

798 

0.049 

ZnS 

2. .342 

40.429 

5.273 

89.272 

2123 

0.081 

ZnSe 

2.454 

32.200 

4.562 

82.687 

1788 

0.080 

ZnTe 

2.637 

29.445 

4.659 

62.430 

1511 

0.071 

CdTe 

2.806 

26.569 

2.722 

38.453 

1371 

0.067 

HgTe 

2.798 

26.396 

2.746 

40.363 

943 

0.056 

■Reference  30. 


AE  =  [3a/(6-6o)V2-f0,(6-8o)V2-(-a(6  +  2y)V2‘ 

+  0(6  +  r,  /2  +  0b^/4  +  0{36  +  k^r)^/S 

8  -f  A-jy ) V4  ]rf  2  \/2Cy V(  3  \/3 )  +  A£ch  . 

(16) 

where  the  chemical  contnbution  is  written  as 

=  4  Ech  <  28 -f  y  )<f  a- V  W ( 28 -f- r  ^ 

-f[t/(AK,)=-G  AEj]  ,  (17) 

where  AE,h  =  E -  IE  A E,  and  AE^  =  t/(  AEj  )^-G  AE,. 
Ech  is  a  chemical  force,  which  when  it  is  positive  tends  to 
push  the  C  atom  away  from  the  impurity  atom  A.  This 
force  arises  from  the  difference  in  the  bond  tensions  in¬ 
duced  between  the  AC  and  BC  bonds  adjacent  to  C  be¬ 
cause  the  neighbsrring  antibinding  states  are  different 
from  those  of  their  respective  host  states.  AE^  is  due  to 
the  difference  in  the  polarities  AE3  alone  and  is  indepen¬ 
dent  of  the  displacement.  Finally,  H  can  be  regarded  as  a 
chemically  induced  force  constant,  which  when  it  is  posi¬ 
tive  tends  to  restrain  the  lattice  from  distortion  and  in¬ 
creases  the  elastic  energy. 

The  equilibrium  requirements  3lAE)/d8  =  0  and 
3(AE)/3y— 0  then  lead  to  the  solution  y  =  Qb,  and  8  is 
given  by 

8  =  (8o -(-8^ )/|  1 -(- [  a(  1  -  2Q) -cE(  17/4  -  Ag) 

-hl6//(l-2(?l/9]/(3a;-(-/3,)t  ,  (18) 

where  the  constant  /.  is  I  -1-  3/.|/4-t-A.;/2,  and 


I 

I 


&o=-4Eeh(l-e/2)/(3£f(3a,-fft)]  ,  (19) 

with  Q  =  2J /K,  J  r^a  +  kP/2  +  %H /9,  and 
K=4a  +  2\^CdA3\^)A-{\+k]/4-yk\/2)P-ym/9  . 


IV.  VALENCE-FORCE-FIELD  MODELS 

In  this  section  we  consider  several  models  based  on  the 
valence  force  field.  These  models  have  been  used  fre¬ 
quently  to  explain  the  impurity  bond  relaxation.’  ’  We 
shall  first  derive  the  explicit  expressions  for  these  models 
and  then  connect  them  with  the  existing  results. 

A.  Model  A:  Third-shell  atoms  and  beyond  are  fixed 
at  their  pure  crystal  positions 

Let  the  bond  lengths  surrounding  the  impurity  again  be 
di=(f'l— 6)  and  let  the  second-jhell  atoms  have  radial 
displacements  of  the  forms  id/v  3),  (y.y.O).  etc.  Beyond 
and  including  the  third  shell,  all  the  other  atoms  are  held 
at  their  pure-crystal  positions.  There  are  nine  different 
contributions  to  the  strain  energy  in  VFF  (see  Appendix 
Bi:  the  a  terms  from  the  four  bonds  surrounding  the  im¬ 
purity,  6a/(8-6o)‘<f’:  the  0  terms  among  the  six  pairs  of 
these  bonds,  /3;(6  — 8o)’if^  the  0  terms  between  the  first- 
and  second-shell  bonds,  2j3l8-(->')’<f’;  the  n  terms  from 
the  second-shell  bonds,  2aib-*-2y  \'d'\  the  0  terms  among 
the  second-shell  bonds,  20b'd'\  the  0  terms  between  the 
second-  and  third-shell  bonds,  j0b'd'  -^0tb  +  2y)'d'-,  the 
a  terms  from  the  third-shell  bonds,  &ay'd':  the  0  terms 
among  the  third-shell  bonds,  40y^d^\  and  the  0  terms  be¬ 
tween  the  third-  and  fourth-shell  bonds,  b0y^d'.  Thus, 
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the  excess  energy  (in  this  case  7  times  the  strain  energy) 
becomes 

A£  =  [  |a,(8-6o)^+T^;(6-6o)^+T^(6  +  y>^ 

+  I  a(  8  +  2 y  -f  ^ /36^  +  7^(  6  +  2y 

+  2ar+jPy^]d^ .  (20) 

The  minimization  of  AE  with  respect  to  8  and  y  leads  Eq. 
(20)  to  y  =  —  8/4,  and 

6=8o/[l+(a+17^/2)/(6a,-)-/3/)]  .  (21) 

We  note  that  there  is  some  ambiguity  in  the  third  con¬ 
tribution  listed  above  for  the  P  terms  between  the  impuri¬ 
ty  and  host  bonds.  The  value  of  0  could  be  chosen  as  one 
of  these  combinations  /3,  ^/,  j{p+0i),  or  other 

proper  combinations.  Because  the  values  of  P  and  are 

comparable  and  /3  values  are  much  smaller  than  a  (see 
Table  1),  the  results  for  8  and  AE  are  not  too  sensitive  to 
the  choice.  There  is  also  some  ambiguity  in  the  values  for 
“undistorted”  crystal.  The  —  d^/3  used 
is  the  simplest  choice.  A  different  choice  will  not  affect 
the  results  for  8  at  all,  but  will  make  AE  slightly  different. 
In  fact,  model  A  was  first  used  by  Martins  and  Zunger.^ 
However,  their  expression  for  8  is  different  from  Eq.  (21) 
because  they  made  different  choices  of  the  two  quantities 
just  mentioned.  Nevertheless,  Sec.  IV  will  show  that  these 
two  expressions  yield  very  similar  results.  These  ambigui¬ 
ties  do  not  occur  in  the  full  theory  in  Sec.  Ill,  where  the 
impurity-host  interactions  are  taken  into  account  natural¬ 
ly  by  the  replacement  energy  A^  (see  Eq.  (10)]. 

B.  Model  B:  Second-shell  atoms  connect 
to  a  fixed  boundary 

This  model  corresponds  to  ;'=0  in  model  A.  So  we 
have 


A£  =  [  ja,(8-8o)^-f-T^/(8-8o)^4- 

(22) 

and 

8  =  8o/ll-)-(a-f-19j3/4)/(3a,-(-|^,)]  .  (23) 

This  expression  will  be  used  to  study  the  effect  of  trunca¬ 
tion. 


C.  Model  C:  Simple  spring  model 

If  all  the  P‘%  in  Eqs.  (22)  and  (23)  are  set  equal  to  zero, 
we  have  the  simple  spring  model  with 

AE  =  [ia,(6-8o)^+iab^]d^  (24) 

and 

8=8o/(l-l-Ta/«/>  •  (25) 

The  spring  model  recently  discussed  by  Shih  et  al}  corre¬ 
sponds  to  Eq.  (25)  with  a=ai,  so  8/80=  7. 

D.  Model  D:  VFF  with  the  continuum  connected 
to  the  second-shell  atoms 

1.  Model  Dl 

In  this  case,  A£  only  contains  the  first  five  contribu¬ 
tions  listed  for  case  A  plus  the  elastic  energy  in  the  con¬ 
tinuum.  However,  the  P  terms  between  the  second-  and 
third-shell  bonds  are  modified  because  atoms  outside  R  in 
the  continuum  now  have  radial  displacements  proportion¬ 
al  to  the  inverse  of  the  square  of  the  radius.  The  result  is 


A£  = 


V2 


Tctjib  —  7^/(8  —  80)*-+-  70[(8-f-2y)^-(-  y£(8-t-y)^-t-  yPb^  -t-7^(38-(-A.jy)^-(-  4)3(8-fA.2y)‘'-(-  Cdy*' 


d' 


(26) 


where  /|  and  A.i  are  the  same  as  the  constants  that  appear  in  Eq.  (16).  The  corresponding  equilibrium  condition  can  be 
shown  to  be 

6--fi„/|  I  -clot  1 -2(?)-t- 19j?/4  ^  1  •c.U,/4^/;/4)^]/(3a, +/?,/2l(  ,  (27) 


where 

Q 


(4cf  +-  2v  2Cd /3v  It  ^ p  *  ^k]P  r  7  Ay?) 
2.  Model  02 


i28) 


jP,(b~bo)^d^+iPb^d^ 

in  Eq.  (16)  become 

{p,(b~btt)‘d'^\pb^d' 


A  coniparison  between  Eu  '26  :init  the  full  pcnuiha. 
tion  theory,  Eq  1I61,  shows  :wo  m.iiiir  Jiffereiices  Eir^t. 
all  the  chemical  terms  ,ire  .ihseni  in  Eq.  26',  Secomily, 
the  terms 


in  Eq.  i26i.  This  difference  in  the  strain  energy  will  mask 
the  true  effeets  of  chemical  forces  if  h  from  Eq.  26  is 
compared  with  EPT  .A  better  w.o  to  study  the  ctiemical 
effects  IS  to  use  the  lollowing  equ.ition 
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A£  =  (ia, (6-5)^+ i/3, (8-8o)^-t-ia(6  +  2y)^+^^(5+r>^+ 7/56^+ 1/5(36  +  3., y)^+i^l(6  +  A.jy)^lrf^+^Cy^d’  . 

(29) 

which  i&  Eq.  (16)  with  all  the  chemical  terms  neglected.  The  corresponding  6  becomes 

6=6o/ll+[a(l-2O)+3ip_0(l  +  iX,+  i3.j)Ql/(3a,+0,)l  .  (30) 


with  Q  still  given  by  Eq.  (28). 


I  - - - — 

where  the  mixing  energy  is  given  by 


E.  Model  E:  Continuum  connected 
to  the  rmt-shell  atoms 

In  this  case,  y  =  —  3t/36/(8t/2)  and  Af  only  includes 
the  first  three  contributions  listed  in  model  A  plus  the 
strain  energy  of  the  continuum: 


A£  = 


ia,(6-8o)^+i^,(6-6o)^ 


1- 


3vl 

8v^ 


y^+icas^ 


(31) 


The  relaxation  parameter  is  ^iven  by 

8=6o/ 


3vl 

1 

/  , 

1  + 

jCd  + 

i 

'  8v^2 

0 

/(3a, + 

(32) 


We  note  that  the  continuum  model  used  to  estimate  the 
bond-length  relaxation  by  Baldereschi  and  Hopfleld’  cor¬ 
responds  to  Eq.  (32)  without  the  /3  terms,  which  yields 
6/6o!«0.4  to  0.3,  rather  than  the  proper  values  around  0.7 
to  0.8. 


A£„  =  ^»noy 

=^6Af(t4AfAA  -^^ab''ab+^bb''bb^~(>^^^^aa  +3'Eaa I 

=  6^rA£,  (35) 

where 

=  ~ +ebb)  •  (36) 

The  mixing  entropy  AS  can  also  be  written  from  a  simple 
generalization  of  the  random  distribution.'^  For  modest 
pressure,  AE  is  the  same  as  the  mixing  enthalpy  AH„ . 

Now  the  |>air  interaction  energies  can  be  approximately 
related  to  the  impurity-substitution  energies  by 


A,(i4  in  BC) as  12(e,4j —  £(„ ) 

(37) 

and 

Ai(A  in  BC)«  l2(E,ig  —egg ) . 

(38) 

Thus,  Ae  of  Eq.  (36)  becomes 

Ae  =  -^(A,(A  in  BC)  +  A,(B  in  AO] 

=  ^(A£(/I  in  BC)  +  A£(B  in  ^O]  . 

(39) 

V.  ALLOY  MIXING  ENTHALPY 

The  impurity-substitution  excess  energies  A£  provide  a 
first  estimate  of  the  mixing  enthalpies  of  pseudobinary  al¬ 
loys.  Most  current  thermodynamics  theories  of  semicon¬ 
ductor  alloys  are  based  on  an  extension  of  the  binary  solu¬ 
tion  model. In  this  model,  the  mixing  Helmholtz  energy 
of  an  alloy  is  defined  as 

=^illoy —  (.^^/IC+Y^Bc)  .  (33) 

where  y  =  1  — x,  and  F^c  ^bc  are  the  respective  free 
energies  of  the  pure  AC  and  BC  compounds  at  the  same 
temperature.  Because  the  C  atoms  occupy  a  sublattice, 
the  nearest  neighbors  of  A  and  B  atoms  in  the  alloy  are 
the  C  atoms.  Thus,  the  pair  potentials  that  enter  the 
binary  solution  theory  are  now  the  second-neighbor  in¬ 
teractions.  Let  N N^s,  and  Ngg  be,  respectively,  the 
numbers  of  the  second-neighbor  AA,  AB,  and  BB  pairs, 
with  corresponding  pair  interaction  energies  b^A’^aB’  and 
egg.  For  tetrahedral  semiconductors,  there  are  a  total  of 
bN  second-neighbor  pairs  for  a  crystal  containing  A'  unit 
cells.  Denote  the  ratios  N^a,  ^Ag,  and  Ngg  to  6N  as 
'‘aa’  ^AB'-f’  and  rgg,  respectively.  Then  those  ratios  are 
related  to  the  alloy  composition  by  r.,^=x— r/2  and 
^BB  =Y  -  r/2.  The  mixing  free  energy  has  two  terms, 

(34) 


Usually,  the  experimental  A//„  is  written  as 

AH„=x(l-x)n,  (40) 

which  is  equivalent  to  assuming  a  random  distribution, 
i.e.,  r=2x(l— x).  Using  this  expression  for  r  and  com¬ 
paring  Eqs.  (40)  and  (35),  we  see  that  the  mixing  enthalpy 
parameter  H  is  given  by 

n  =  2lA£(A  in  BC)  +  A£(B  in /lO]  .  (41) 

This  connection  provides  a  further  check  of  the  theory. 

VI.  NUMERICAL  RESULTS  AND  DISCUSSION 
A.  Chemical  terms 

Table  11  lists  6(,=  l-d,/d,  d=  \  —d^/d,  the  excess  en¬ 
ergy  (per  bond)  A£  for  the  full  theory  and  its  correspond¬ 
ing  VFF  model  D2,  and  the  terms  derived  from  the 
metallization  energies.  6^  [Eq.  (19)],  F^h.  H,  ^Ep,  and 
A£ch  (Eq.  (17)).  The  appropriate  derivatives  /,  g,  h,. . . 
(see  Eq.  (9)]  are  computed  using  the  atomic-term  values 
that  we  have  generated  from  impurity -level’’  and 
structural  studies.- 

For  substitutions  involving  the  cation  pair  (Ga,A)),  F^^ 
has  the  same  sign  as  6(,,  which  means  that  prevents 
relaxation  and  thus  tends  to  increase  the  strain  energy. 
The  chemical  forces  H  are  also  significant.  As  a  result. 
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TABLE  II.  Comparison  between  the  full  theory  and  the  corresponding  VFF  model  D2  to  study  the  effects  of  chemical  terms.  All 
A£*s  are  in  units  of  kcal/moleband. 


Model  02 


Full  theory 


Impurity 

Host 

fh 

6 

\E 

61. 

6 

fth  <10''®N) 

H  (N/m) 

Ga 

AlP 

0.003 

0.002 

0.001 

-0,001 

0.001 

0.150 

4,581 

-0.016 

-0.016 

-0.013 

At 

GaP 

-0.003 

-0.002 

0.001 

0.001 

-0.001 

-0.150 

4.581 

-0.016 

-0.016 

-0.013 

Ga 

AlAs 

0.001 

0.001 

0.000 

-0.002 

-0.001 

0.243 

5.733 

-0,020 

-0.021 

-0.018 

A1 

GaAs 

-0.001 

-0.001 

0.000 

0,002 

0.001 

-0.243 

5.733 

-0.020 

-0.021 

-0.018 

Ga 

AlSb 

0.006 

0.004 

0.005 

-0,004 

0.001 

0.389 

5,632 

-0.054 

-0.053 

-0.039 

Al 

GaSb 

-0.006 

-0.004 

0.005 

0.004 

-0.001 

-0.389 

5.632 

-0.054 

-  0.053 

-0.039 

In 

GaP 

-0,077 

-0.052 

0.959 

-0.006 

-0.054 

0  699 

3.778 

-0.188 

-0.219 

0.742 

Ga 

InP 

0.071 

0,056 

0.734 

0.005 

0.057 

-0.699 

3.778 

-0.188 

-0.206 

0.530 

In 

GaAs 

-0.071 

-0.048 

0.752 

-0.009 

-0.050 

0.804 

4.778 

-0.257 

-0.283 

0.472 

Ga 

InAs 

0.066 

0.052 

0.592 

0,007 

0.054 

-0.804 

4.778 

-0.257 

-0.265 

0.330 

In 

GaSb 

-0.062 

-0.043 

0.554 

-0.004 

-0.042 

0.352 

5.201 

-0.363 

-0.308 

0.247 

Ga 

InSb 

0.059 

0.046 

0.445 

0.004 

0.044 

-0.352 

5.201 

-0.363 

-0.287 

0.160 

In 

AlP 

-0.074 

-0.053 

0.761 

-0.007 

-0.056 

0.769 

3  506 

-0.035 

-0.087 

0.679 

Al 

InP 

0.068 

0.053 

0.674 

0.006 

0.056 

-0.769 

3.506 

-0.035 

-0.083 

0.596 

In 

AlAs 

-0.070 

-0.048 

0.705 

-0.010 

-O.OSl 

0.942 

4.43? 

-0.048 

-O.lll 

0.602 

Al 

InAs 

0.065 

0,052 

0.576 

0.008 

0.054 

-0.942 

4.437 

-0,048 

-0.099 

0.485 

In 

AlSb 

-0.056 

-0.0.39 

0.440 

-0.008 

-0.041 

0.689 

4.979 

-0.061 

-0.073 

0.369 

Al 

InSb 

0.053 

0  042 

0.,368 

0.007 

0.044 

-0.689 

4.979 

-0.061 

-0.061 

0.310 

Cd 

ZnTf 

-0.064 

-0.048 

0.432 

-  0.003 

-0.050 

0.202 

--  0.484 

-0.005 

-0.064 

0.373 

7.n 

CdTe 

0.060 

0.050 

0.314 

0.002 

-0.053 

-0.202 

-0.484 

-0.005 

-0.072 

0  247 

Hg 

CdTe 

0.003 

0.002 

0.001 

0.004 

0.005 

-0.278 

-0.753 

-0.018 

-0.026 

-0018 

Cd 

HgTe 

-  0,003 

-0,002 

0.001 

-0.004 

-0.005 

0.278 

-0.753 

-0.018 

-0.026 

-0.018 

Hg 

ZnTe 

-0  061 

-0.045 

0.392 

-0,001 

0.046 

0.075 

0.002 

0,052 

0.037 

0.429 

Zn 

HgTe 

0.058 

0.048 

0.286 

0.001 

0.049 

-0.075 

0.002 

0.052 

0,035 

0.322 

As 

AlP 

-0.035 

-0.026 

0.179 

0.001 

-0.025 

-0.085 

0.717 

-0.005 

0.008 

0.187 

P 

AlAs 

0.0.34 

0.025 

0.185 

-0.001 

0.025 

0.085 

0.717 

-0.005 

0.008 

0.194 

As 

GaP 

-0.037 

-0.025 

0.226 

0,002 

-0.024 

-0.181 

1.078 

-0.011 

0.012 

0.240 

P 

GaAs 

0.036 

0.027 

0.21 1 

-0.001 

0.025 

0.I8I 

1.078 

-0  011 

0.014 

0.228 

As 

InP 

-0  032 

-0  023 

0  136 

0.001 

-  0.022 

-0.057 

0.919 

-  0.003 

0,008 

0.144 

P 

InAs 

0.031 

0.024 

0  128 

0,001 

0.024 

0.057 

0.919 

-0.003 

0  009 

0  1.38 

Sh 

AlAs 

-0.084 

-  0,058 

1.024 

0.008 

-0.051 

-0815 

0,644 

-0.180 

0.(X)2 

l.ObO 

As 

AlSb 

0.077 

0.059 

0.919 

-0.007 

0.053 

0  815 

0.644 

-0.180 

0027 

0.984 

Sb 

GaAs 

-  0,078 

0.052 

0.908 

0.018 

-0.040 

-  1  599 

0,927 

0  363 

0.106 

0  929 

As 

GaSb 

0,073 

0.055 

0.823 

-0,014 

0.044 

1.599 

0.927 

-  0.363 

0061 

0  904 

Sb 

InAs 

-0,070 

-0.051 

0.603 

0,010 

-0.042 

-0.824 

0.855 

0.171 

0.009 

0.645 

As 

InSs 

0.065 

0,051 

0.551 

-0.009 

0.044 

0,824 

0  855 

-0  171 

0.CX)8 

0.613 

Sh 

AlP 

-  0.122 

-0,085 

2.007 

0.010 

-0,077 

-0  944 

0  645 

0.241 

0.074 

2,127 

P 

AlSb 

0.100 

0.085 

1.855 

-  0.008 

0,078 

0  944 

0  645 

-0.241 

0.123 

2,0,30 

Sb 

Gap 

-0  no 

-0.075 

2.1.32 

0.021 

-0.061 

-  1  868 

0  930 

0.50S 

0  046 

2  244 

P 

<  iaSh 

0.106 

0,083 

1.806 

-0.015 

0.070 

1.868 

0  9  30 

0  505 

0,09.3 

2.1184 

Sb 

InP 

0  uu 

0,072 

1  383 

0,01 1 

0.063 

t>9;2 

0  854 

0,214 

0.059 

1  501 

P 

InSb 

0  004 

0,077 

1.191 

0.(X)8 

0  069 

0922 

0  854 

0.214 

II  123 

1  37'i 

Se 

ZnS 

0,048 

-0  036 

0.231 

OOOl 

0.036 

0  07' 

0  615 

(1  (Xl.' 

(1  (XX) 

0  2.1 1 

S 

ZnSe 

0  046 

0.037 

0.221 

0  001 

00.37 

0.0' 7 

0,645 

0,(X13 

'1 IXI] 

0  222 

Tc 

ZnSe 

-0,075 

0  056 

0.550 

0.000 

-  0.056 

0  028 

0  6.35 

II  (XI8 

U  1124 

0  s74 

Se 

ZnTe 

0.060 

0.054 

0  532 

0.«X) 

0,054 

-0.028 

0.6.35 

0  (XIK 

0.025 

0  557 

3 
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all  six  cases  involving  this  pair  have  nearly  equal  d  i  and 
d;,  t  e.,  the  small  bond-length  differences  are  made  even 
smaller  in  the  alloy.  The  excess  energies  all  become  nega¬ 
tive,  mainly  because  A£p  is  negative.  For  the  systems  in¬ 
volving  the  (Ga.ln)  and  (In,Al)  pairs,  has  the  opposite 
sign  from  &o,  so  6o  and  6o  have  the  same  sign.  The  chemi¬ 
cal  force  favors  bond  distortion.  However,  because  H  is 
positive  and  introduces  an  increase  in  the  denominator  of 
Eq.  (18),  most  of  the  effect  of  bp  is  cancelled.  For  cases 
involving  (Ga,In),  the  polarity  contributions  A£p  are  all 
negative.  The  fd,  Ad  term  is  negative,  but  //(Ad)^  is  pos¬ 
itive,  so  they  cancel  to  a  certain  degree  and  leave  A£ 
lowered  primarily  because  of  A£p.  While  AEp  is  still  neg¬ 
ative  for  the  (ln,Al)  substitutions,  its  magnitude  is  reduced 
considerably.  The  other  chemical  energies 
fch  Ad -)-//( Ad can  be  as  large  as  AEp,  but  the  overall 
reductions  of  A£  are  only  a  fraction  of  those  for  the 
(Ga,In)  cases.  For  the  several  Il-VI  systems  studied,  both 
fch  and  H  are  small  and  the  net  changes  in  b  have  the 
same  sign  as  bp.  However,  because  bp  is  small  in  the 
(Cd,Hg)  substitutions,  £<.h  actually  causes  a  reversal  in 
which  the  short  bond  length  becomes  shorter  and  the 
longer  one  becomes  longer.  This  is  the  only  exceptional 
case  of  this  type  found  for  all  the  systems  studied.  The 
change  of  A£  due  to  chemical  terms  in  the  (Hg,Zn)  substi¬ 
tution  IS  also  peculiar — it  increases  mainly  because  AEp  is 
positive. 

Next,  we  examine  the  anion  substitutions.  For  the 
groups  involving  the  (P,As)  pair,  the  chemical  shifts  are 
all  small,  but  the  trend  is  less  toward  relaxation  and  larger 
A£.  The  groups  involving  (As,Sb)  and  (P,Sb)  pairs  behave 
very  similarly:  are  significant  and  are  opposing  relax¬ 

ations,  i.e.,  bi  and  bp  have  opposite  signs.  At  the  same 
time,  the  H  values  are  several  times  smaller  than  those  for 
the  corresponding  Ill-V  cation  substitution  case.  Thus, 
most  of  bp  translates  into  a  real  reduction  of  the  ratio 
b/6(i  and  consequently  introduces  extra  strain  energy.  Al¬ 
though  the  AEp  energies  are  significant  and  negative, 
F,.),  Ad  are  positive  and  the  net  A£(.h  can  be  either  positive 
or  negative.  However,  the  induced-strain  energy  due  to 
reduction  of  the  b/b„  makes  all  A£  positive  for  these  two 
groups  of  systems.  For  11-VI  systems,  all  the  chemical  ef¬ 
fects  again  are  small,  but  the  net  chemical  changes  on  AE 
are  slightly  positive. 

The  above  discussion  can  be  summanzed  in  Fig.  2, 
where  the  excess  energies  AE  calculated  from  the  full  per¬ 
turbation  theory  and  MixJel  D2  are  plotted  normalized  to 
the  results  of  the  simple  spring  mvxlel  of  Shih,  Spicer, 
Harrison,  and  Sher  <SSHS)  iRef,  8),  i.e..  Eq.  f24i  with 
a/  =a,  so  A£/(  ,  rrt/‘!  =  bn.  The  calculated  AE  nses  fas¬ 
ter  for  b() '0  than  for  bp  10,  mainly  because  a/a/  =  I.  In 
fact,  if  the  relation’  ’®  a/af  =  \d[/d)^  with  5  of  order  of 
3  to  5  is  used  in  Eq.  >24),  we  obtain  a  percentage  correc¬ 
tion  of  b’bp/4  to  the  SSHS  results,  which  explains  the 
skewed  behavior  of  the  curve.  It  is  also  clear  that  AE 
rises  faster  than  b,',  for  larger  b,,.  However,  the  zeroth- 
order  theory  of  SSHS  is  clearly  an  excellent  representation 
of  the  global  features  of  AE.  The  results  from  model  D2 
are  closer  to  the  parabolic  form  than  those  from  FPT 
The  figure  clearly  shows  the  general  trends,  the  chemical 
terms  cause  negative  shifts  in  AE  for  cation  substitutions 


FIG.  2.  The  excess  energies  A£  over  jad^  calculated  from 
the  full  perturbation  theory  (FPT)  and  its  corresponding 
va'ence-force-field  model  D2.  The  solid  curve  corresponds  to 
Eq.  (24)  with  a  =  a/. 


and  positive  shifts  for  anion  impurities.  It  is  also  clear 
that  the  chemical  shifts  can  be  very  large.  These  effects 
will  have  important  consequences  on  the  alloy  mixing 
enthalpies  to  be  discussed  later. 


B.  Impunty  bond  length 

Tabic  III  lists  the  impunty  bond  lengths  d,  calculated 
from  different  models,  while  a  comparison  of  theory  and 
the  available  experimental  data'  is  presented  in  Table 
IV.  The  actual  size  of  changes  in  d  |  induced  by  the 
chemical  terms  can  be  seen  by  comparing  model  D2  with 
the  full  theory  Except  for  the  systems  involving  the  sub¬ 
stitution  of  iGa.Asi  and  'P.Sbi  pairs  'where  changes  range 
from  0.01  to  0.03  A',  all  the  chemically  induced  changes 
are  less  than  0.01  A  Companson  among  mrxlels  A.  B, 
and  C  shows  that,  while  extending  the  boundary  helps  the 
relaxation  'compare  model  B  to  Ai,  i.e.,  b^b,  is  closer  to 
I.  the  inclusion  of  the  Nind-bendmg  forces  'the  P  terms' 
'compare  miKlel  B  with  f  prevents  it  The  simple  spring 
nuxiel  'model  C'l,  which  contains  neither  o(  these  terms, 
evidently  represents  a  delicate  cancellation  of  these  effects 
and  predicts  results  clos^-  to  those  of  the  lull  perturbation 
theory  and  expenment  .Although  the  d]  values  of  model 
C  are  often  very  close  to  those  of  mixJel  A.  there  are  cases 
(e  g  ,  Ga'P.Sb']  tn  which  mixlcl  C  can  differ  from  model 
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TABLE  III.  Calfulaled  impurity  local  bond  lengths  (in  A)  from  the  full  theory  and  several  valence 
force  mixlels  discussed  in  Sec.  Ill,  and  their  comparison  with  the  values  calculated  by  Martins  and 

Zuriger  (Ref.  3).  _ _ _ _  _  _ _ 

rmpurity  ^  A  B  C  Dl  D2  E  FPT  MZ 


Ga 

AIP 

2.362 

2.363 

AI 

GaP 

2.365 

2.364 

Ca 

AlAs 

2.449 

2.449 

Al 

GaAs 

2.450 

2.450 

Ga 

AlSb 

2.645 

2.646 

Al 

GaSb 

2,651 

2.650 

In 

GaP 

2.477 

2.462 

Ga 

InP 

2.406 

2.421 

In 

CaAs 

2.559 

2.544 

Ga 

InAs 

2.492 

2.506 

In 

GaSb 

2.747 

2.734 

Ga 

InSb 

2.6g3 

2.697 

In 

AIP 

2.4«7 

2.472 

Al 

InP 

2.412 

2.427 

In 

AlAs 

2.561 

2.546 

Al 

InAs 

2.493 

2.506 

In 

AlSb 

2.754 

2.741 

Al 

InSb 

2.693 

2.705 

Cd 

ZnTe 

2.756 

2.740 

Zn 

CdTe 

2.673 

2.688 

Hg 

CdTe 

2,800 

2.801 

Cd 

HgTe 

2.804 

2,803 

Hg 

ZnTe 

2.750 

2.735 

Zn 

HgTe 

2.671 

2,685 

C 

Dl 

D2 

E 

FPT 

2.362 

2.362 

2.362 

2.364 

2.364 

2.365 

2  365 

2.365 

2,365 

2.363 

2.449 

2.449 

2.449 

2.450 

2.452 

2.450 

2.450 

2.450 

2.450 

2.447 

2.644 

2.645 

2.644 

2.649 

2.653 

2.652 

2.652 

2.652 

2.650 

2.644 

2.492 

2.479 

2.483 

2.435 

2.487 

2.402 

2.395 

2.399 

2.409 

2.396 

2.573 

2.561 

2.565 

2.518 

2.570 

2.486 

2.482 

2.485 

2.496 

2.481 

2.760 

2.749 

2.754 

2.710 

2.750 

2.678 

2.673 

2.676 

2.686 

2.680 

2.494 

2.490 

2.493 

2,447 

2.498 

2.408 

2.401 

2.405 

2.415 

2.400 

2.572 

2.563 

2.568 

2.523 

2.575 

2.487 

2.483 

2.487 

2.497 

2.480 

2.763 

2.756 

2.760 

2.721 

2.765 

2.689 

2.685 

2.687 

2.696 

2.683 

2.760 

2.760 

2.763 

2.720 

2.770 

2.676 

2,660 

2.665 

2  671 

2.658 

2.800 

2.800 

2.799 

2.801 

2,790 

2.804 

2.805 

2.805 

2.804 

2.813 

2.754 

2.753 

2.757 

2.715 

2.758 

2.674 

2.659 

2.664 

2.671 

2.662 

•w' 

As 

AIP 

2.425 

2.418 

2.429 

2.427 

2-428 

2  406 

2.427 

2.422 

P 

AlAs 

2.392 

2.399 

2.387 

2.387 

2.389 

1.394 

2.390 

2  395 

■1 

As 

GaP 

2.417 

2.409 

2.424 

2.417 

2,420 

2  396 

2.416 

2  414 

P 

GaAs 

2.386 

2.393 

2.380 

2.382 

2.383 

2.389 

2  386 

2,387 

r_- 

As 

InP 

2.596 

2.589 

2.599 

2.598 

2.600 

2.579 

2.598 

2.595 

P 

InAs 

2.561 

2.568 

2.558 

2.557 

2.558 

2.563 

2.560 

2  562 

f.-* 

Sb 

AlAs 

2.584 

2.566 

2  .597 

2.587 

2.592 

2  539 

2  577 

2,574 

As 

AlSb 

2.506 

2.522 

2.496 

2.495 

2,498 

2  511 

2  514 

2  510 

Sb 

GaAs 

2.569 

2.553 

2.584 

2.571 

2,576 

2,524 

2  546 

2  564 

As 

GaSb 

2.501 

2.516 

2.489 

2.492 

2  495 

2  508 

2  525 

2  '05 

Sb 

InAs 

2.747 

2.730 

2.754 

2.750 

2.754 

2.705 

2  '.33 

2  739 

As 

InSb 

2.669 

2.683 

2.663 

2.658 

2  662 

2  672 

2  68  3 

2  66' 

L: 

Sb 

AIP 

2.555 

2.529 

2  569 

2  561 

2  568 

2  488 

2  550 

2  '42 

P 

AlSb 

2.440 

2.462 

2  426 

2  425 

2  4. to 

2  44' 

2  448 

2  444 

* 

Sb 

GaP 

2.526 

2  503 

2  551 

2.529 

2  537 

2  46! 

2  504 

2  '|9 

P 

GaSb 

2.431 

2.451 

2.414 

2  418 

2  422 

2  440 

2  4.S4 

2  43t> 

Sb 

InP 

2.712 

2.687 

2.720 

2.719 

2  725 

2  654 

2  '02 

2  'IKI 

C%’ 

P 

InSb 

2.599 

2.619 

2. .'91 

2-585 

2  590 

2  604 

2.611 

2  59' 

fe 

Se 

ZnS 

2  420 

2.409 

2.421 

2.424 

2.426 

2  39(1 

2  426 

2  420 

S 

ZnSe 

2.367 

2.376 

2.365 

2  360 

2  364 

2  .370 

2  36.3 

2  36’ 

Te 

ZnSc 

2.586 

2  569 

2  588 

2.589 

2  592 

2  .540 

2  .‘9| 

2  's4 

• 

Sc 

ZnTe 

2.. 501 

2.5  r 

2.497 

2.490 

2  494 

2  504 

2  49“^ 

2  502 

Tc 

ZnS 

2.543 

2.513 

2.544 

2.552 

2.558 

2  478 

'  s  s-* 

2  '  39 

S 

ZnTe 

2.406 

2  429 

2.4(X) 

2. '90 

2  .396 

2  410 

2  .'-''n 

2  40' 

It-i 
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A  by  0.025  A.  Model  A  produces  about  the  same  di 
values  as  model ^Dl,  where  the  maximum  difference  in  dj 
is  only  0.015  A.  Martins  and  Zunger'  used  the  same 
mode)  as  model  A;  however,  their  analytic  expression  for 
6  is  different  from  that  given  by  Eq.  (21).  Nevertheless, 
numerical  results  indicate  that  these  two  calculations 
agree  to  0.01  A.  The  slightly  different  forms  of  strain  en¬ 
ergies  used  in  models  D2  and  DI  only  introduce  a  small 
change  in  <f|,  with  the  largest  difference  being  less  than 
0.01  A.  The  first-shell  continuum  model  (model  E)  allows 
too  little  relaxation,  so  while  the  other  models  produce  a 
ratio  6/6o  ranging  from  0.6  to  0.8,  model  E  only  ranges 
from  0.4  to  0.6.  The  reason  that  the  fixed  boundary  in 
model  A  works  is  that  the  effective  shear  coefficient  C 
(see  Table  1)  characterizing  the  strain  energy  in  the  elastic 
continuum  is  large.  However,  model  B  is  too  rigid  and 
does  not  provide  enough  buffer  between  the  impurity 
bond  and  the  fixed  boundary. 

The  comparison  of  the  theoretical  results  with  the 
available  experimental  data  in  Table  IV  indicates  that 
models  B  and  E  are  the  least  accurate.  Models  A,  DI, 
that  of  Martins  and  Zunger  (MZ),  and  the  full  theory  are 
comparabje  in  that  all  have  an  average  absolute  deviation 
of  0.012  A,  which  is  about  the  experimental  uncertainty  in 
extended  x-ray-absorption  fine  structure  (EXAFS).  The 
agreement  between  theory  and  experiment,  however,  is  not 
uniform.  The  most  surprising  result  in  Table  IV  is  that 
the  simple  spring  model  (model  C)  and  its  cruder  version 
used  by  Shih  et  al}  {a=ai  so  6/6o=0.75,  labeled  as 
SSHS)  have  the  smallest  variance  in  d|,  about  0.006  A. 
We  know  this  does  not  imply  that  the  simple  spring 
model  represents  the  real  picture  of  bond-length  relaxa¬ 
tion.  For  example,  if  we  let  all  the  shear  coefficients  be  0, 
i.e„  0=C  =  O  in  our  model,  then  as  the  range  of  the  boun¬ 
dary  is  gradually  extended,  the  local  bond  length  will 
eventually  relax  to  the  impurity  bond  length  d|  =  d/,  or 
B  =  bfy  This  can  be  seen  in  model  A  from  Eq.  (211,  where 
6  reduces  to  6o/(l -f  a/bOf ),  rather  than  6o/(  1 -t-a/3a; ) 
as  predicted  by  model  C,  and  in  model  D  from  Eqs.  (27) 
and  (30),  5  becomes  6o>  'he  continuum  is  taken  to  be 
shear'iess.  Considering  that  various  effects  are  included 
that  may  mask  the  absolute  accuracy  of  d)  predictions 
(e.g.,  while  low-tempierature  bond  lengths  are  used  in  the 
calculation,  the  room-temperature  values  of  clastic  con¬ 
stants  are  adopted),  the  agreement  of  various  models  with 
experiments  in  Taole  IV  should  be  regarded  as  excellent. 
There  are,  however,  many  other  impurity  systems  in 
which  the  simple-spring-model  predictions  differ  consid¬ 
erably  from  other  models,  as  is  shown  in  Fig.  3,  where  6  is 
plotted  against  6o  for  the  full  theory.  Those  points  that 
deviate  significantly  from  the  0.75-slope  line  are  the  sys¬ 
tems  with  (As,Sb)  or  (P,Sb>  substitutions.  Additional  EX¬ 
AFS  measurements  on  these  systems  are  needed  to  test 
these  predictions. 


C.  Mixing  enthalpies 

Table  V  lists  the  mixing  enthalpy  parameters  SI  (in 
kcal/mole)  for  a  number  of  alloys  estimated  from  Eq.  (41 1 
for  all  the  models  considered,  along  with  other  theoreti¬ 
cal''*’' and  expenmental  values."’  ''*  As  already  dis- 
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TABLE  V.  Mixing  enthalpy  parameter  fl  tin  kcal/molcl  estimated  from  the  full  perturbation  theory  anj  several  valence  force 
models  discussed  in  Sec.  Ill,  and  companson  with  expenments  and  other  theories. 


A 

B 

C 

Dl 

D2 

E 

FPT 

MZ‘ 

DL" 

FM‘ 

VV' 

Expt.' 

(Ga,Al)P 

0.00 

0.01 

0.00 

0.01 

0.01 

0.01 

-0.05 

(Ga»Al)As 

0.00 

000 

0.00 

0.00 

0.00 

0.00 

-0.07 

002 

0.02 

0.03 

0.11 

0.0 

(Ga.ADSb 

0.02 

0.03 

0.02 

0.02 

0.02 

0.03 

-0.15 

0.02 

0.02 

0.03 

0.0 

(Ga.lnlP 

3.76 

4.79 

3.0 

3.29 

3.39 

5.24 

2.54 

4.56 

.V63 

2.94 

3.25,3.5 

IGa,ln)As 

2.97 

3.76 

2.36 

2.61 

2.69 

4.14 

1.60 

2.49 

2.81 

2.42 

1.25 

1.65,  2.0,  3.0 

(Ga,In)Sb 

2.22 

2.83 

1.77 

1.95 

2.00 

3.09 

0.81 

2.53 

1.85 

1.83 

1.47,  1.9 

(In.Al)P 

3.24 

4.22 

2.77 

2.78 

2.87 

4  60 

2.55 

(In.ADAs 

2.86 

3.65 

2.32 

2.49 

2.56 

3.93 

2.17 

3.60 

2.81 

2.37 

2.5 

(In.Al)Sb 

1.81 

2.33 

1.49 

1.57 

1.61 

2  50 

1.36 

2.06 

1.46 

1.45 

0.6 

(Cd,Zn)Te 

1.80 

2.43 

1.73 

1.43 

1.49 

2.45 

1.24 

2.12 

1.97 

1.63 

1.34 

(Hg,Cd)Te 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

-0.07 

0.7,  1.4 

(Hg,2n)Te 

1.63 

2.20 

1.56 

1.30 

1.36 

2.23 

1.50 

1.91 

1.81 

1.48 

3.0 

Al(P.As) 

0.81 

1.03 

0.65 

0.71 

0.73 

1.14 

0.76 

Ga(P.As) 

0.95 

1.18 

0.70 

0.86 

0.87 

1  32 

0.94 

1.15 

0.98 

0.66 

0.12 

o 

b 

In(P,As) 

0.60 

0.78 

0.52 

0.51 

0.53 

0.84 

0.57 

0.72 

0.58 

0.52 

0.4 

Al(As.Sb) 

4.31 

5.45 

3.38 

3.80 

3.88 

592 

409 

Ga(As.Sb) 

3.77 

4.69 

2.81 

3.40 

3.46 

5.22 

3.67 

4.58 

3.35 

2.76 

4.0,  4.5 

In(As.Sb) 

2.61 

3.39 

2.23 

2.24 

2.31 

3.67 

2.52 

2.89 

2.29 

2.17 

6.65 

2.25,  2.9 

AI(P.Sbl 

8.60 

10.99 

6.99 

7.54 

7.73 

12.00 

8.32 

GalP.Sb) 

8.54 

10.61 

6.36 

7.72 

7.88 

11.66 

8.66 

In(P,Sb) 

5.87 

7.64 

5.08 

4.99 

5.15 

8.04 

5.76 

Zn(S,Se) 

1.04 

1.39 

0.98 

085 

0.90 

1.49 

0.90 

ZntS.Sc) 

1.04 

1.39 

0.98 

0.85 

0.90 

1.49 

0.90 

ZntSe.Tel 

2.47 

3.27 

2.23 

209 

2.16 

3.63 

2.26 

2  91 

3.11 

2.12 

3.12 

1.55 

ZntS.Te) 

7.02 

9  34 

6.45 

5.80 

6.02 

9.72 

6  20 

•Reference  3, 

column  A  of  Table  III. 

•Tleference  16. 
'Reference  18. 
‘'Reference  17. 


•References  1 6  and  19. 


cassed,  the  chemical  terms  reduce  the  excess  energies  in 
the  cation  impurities  and  increase  them  for  anion  impuri¬ 
ties.  The  corresponding  changes  in  II  are  the  differences 
between  the  FPT  and  model  D2.  We  note  that  the  reduc¬ 
tions  of  n  for  the  (Ga.In'  alloys  are  very  large  (>l 
kcal/molei  and  also  significant  for  iln.Aj)  alloys.  Howev¬ 
er.  the  increases  in  II  for  the  anion  substitutional  alloys 
are  not  as  large.  Also,  the  II  differences  between  models 
D1  and  D2  are  less  than  10<^f  Mcxiel  A  prtxluces  H 
values  about  20^^  larger  than  nnxiel  Ol,  nnxiel  B  in  turn 
is  20f5-  higher  than  model  A.  and  misJei  E  is  lO'T-  higher 
than  model  B.  The  12  values  in  ihe  simple  spring  model 
(model  Cl  are  seen  t.,  .-e  about  the  same  .is  nuxlcl  Dl.  .’1 
though  the  diffcnnccs  amc'iig  system  can  be  positive  or 
negative  Although  M/.  used  the  same  strain  model  as 
model  A,  their  11  valui-s  do  noi  agree  wiih  our  model  A 
values  because  their  way  of  estimating  II  is  different.  In 


fact,  the  values  of  MZ  are  closer  to  model  B  than  to  A. 

To  distinguish  the  quantitative  nature  of  different 
theoretical  models,  we  note  that  there  are  also  important 
factors  that  may  mask  the  comparison  between  theory 
and  experiment  for  11.  One  important  factor  is  that  the 
mixing  enthalpies  extracted  from  phase-diagram  analysis 
are  sensitive  to  sample  and  expenmental  conditions. 
These  contain  contnbutions  from  vanous  nonideal 
structures  such  as  vacancies,  impurities,  dislocations, 
grain  boundanes,  and  surface  cundiiions,  in  addition  to 
the  ideal  -i//,,  considered  here  for  solid  solutions.  Thus, 
our  theoretical  should  represent  a  lower  bound 

Another  complication  comes  from  the  version  of  the 
theiiry  ot  solid  solution  adopted.  1  he  thesiry  used  for 
analysis  so  far  assumes  a  regular  solid  solution  with 
second -neighNir  pair  interactions  as  was  outlined  in  Sec 
V.  Recent  e.xperimental*'  and  theoretical’’  studies 
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FIG.  4  Pint  ol  /  /  •  ;is  a  fuin.'lion  K  for  /,  oblamfil  tri'm  a  M’  l  ami  h  itiv  i'\pt  tirm.’mal  !!  in  I  ahU  '  In  a',  llio  cir¬ 

cles  are  for  anion  alloss,  the  squares  for  cation  subsiiiiinon.  I  he  solid  lines  in  IvmIi  pans  coircspoml  lo  the  siinple  tlicors  discussed  in 
Appcndi.s  C  The  dashed  lines  al  I,  I  0  separate  the  miscible  from  immiscible  ermips 


FIG.  3.  Calculated  bond-relaxation  parameter  8  from  FPT  and  D2  as  a  function  of  bo-  The  6=780  curve  corresponds  to  the 
theory  of  Shih,  Spicer.  Harrison,  and  Sher  (Ref.  8). 

have  suggested  the  possibility  of  compositional  clustering  composition  variation  of  the  alloy  hardness^^  and  from 
or  long-range  correlations  in  alloys.  Extending  the  theory  the  optical-phonon  frequencies*^  that  the  shear  coeffi-  e's 

to  include  such  effects  will  alter  the  simple  results  for  cients  of  alloys  increase  near  the  center  of  the  composition  ^  E 

AW„  from  Eq.  (411.  Moreover,  there  is  evidence  from  the  range.  This  will  cause  the  effective  continuum  shear  coef- 
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TABLE  VI.  Comparison  of  the  crilical  temperature  T,  of  mixing  and  melting  temperatures  of  the 
constituents  7’,  and  Ti.  in  the  order  of  their  appearance  in  the  parentheses.  Also  tabulated  are  the 
averaged  absolute  values  of  by  ;  and  the  ratio  ;  6,,  ■  /fc„. 


System 

:  hoi 
(%) 

F, 

IK) 

Tc/T, 

TJT, 

'  /bm 

{A].Ga)Sb 

0.6 

0 

0 

0 

0.067 

(Al,Ga)As 

0.1 

0 

0 

0 

0.009 

(AI.Ga)P 

0.3 

0 

0 

0 

0.029 

(Al.In)Sb 

5.5 

342 

0.25 

0.42 

0.679 

(Ga,In)Sb 

5.8 

204 

0.21 

0.25 

0.716 

IAl,In)As 

6.8 

547 

0.29 

0.45 

0.687 

(Ga,In)As 

6.9 

403 

0.23 

0.34 

0.697 

(Al,In)P 

7.1 

642 

0.36 

0.48 

0.732 

(Ga,In)P 

7.4 

639 

0.43 

0.48 

0.763 

In(P,As) 

3.2 

144 

0.11 

0.12 

0.330 

GaiP.As) 

3.7 

236 

0.14 

0.15 

0.352 

Ga(P,As) 

3.5 

191 

0.10 

0.11 

0.307 

In(As,Sb) 

6.8 

635 

0.52 

0.79 

0,840 

Ga(As,Sbl 

7.6 

924 

0.53 

0.94 

0.844 

AllAs,Sb) 

8.1 

1030 

0.56 

0.78 

0.810 

In(P.Sb) 

9.9 

1450 

1.08 

1.82 

1.222 

Ga(P,Sb) 

II. 3 

2180 

1.25 

2.21 

1.256 

AKP.Sbl 

11.6 

2095 

1.19 

1.58 

1.116 

(Cd,HglTe 

0.3 

0 

0 

0 

0.033 

(Zn,Hg)Te 

6.0 

377 

0.25 

0.40 

0.659 

tZn,Cd)Te 

6.2 

312 

0.21 

0.23 

0.564 

Zn(S,Se) 

4.7 

226 

0.11 

0.13 

0..362 

ZniSe.Te) 

7.2 

569 

0.32 

0.38 

0.615 

ZniS.Te) 

11.9 

1560 

0.74 

1.03 

1,017 

ficient  C  in  Eq.  1 14)  to  be  composition  dependent,  which  tendency  toward  compound  formation.  Secondly,  the 

will  cause  fl  to  increa.se.  Despite  these  uncertainties,  use-  FPT  tends  to  predict  smaller  fl  values  than  observed  ex- 

ful  comparisons  across  the  board  in  Table  V  can  still  be  perimentaJly,  which  should  be  expected  according  to  our 

discussion.  To  the  extent  that  the  FPT  predicts  the 
Based  on  the  above  considerations,  we  can  conclude  correct  values  for  an  ideal  solution,  the  difference 

that  models  B,  E,  and  MZ  predict  fl  values  that  are  too  may  be  attributed  to  imperfect  conditions 

high.  We  should  emphasize  that  all  the  fl  numbers  for  and  deviations  from  the  ideal  solution  theory  , 
models  from  A  through  MZ  are  directly  calculated  Finally,  the  calculated  fl  values  in  Table  V  provide 
without  any  adjustable  parameters.  The  fact  that  mtxJels  some  guidance  in  separating  the  completely  miscible  al- 

A,  C,  and  D1,D2,  and  the  FPT  agree  with  the  experiment  loys  from  immiscible  ones.’*'  ’  In  a  true  random  alloy, 

'  as  well  as  or  even  slightly  better  than  the  one-parameter  the  cntcna'"  for  alloy  mixing  at  a  tem[ierature  T  is  that 

theories,  the  delta-lattice  (DL)  model'*’  and  the  model  of  T  '  where  the  critical  temperature  T,.  is  given  by 

Fedder  and  Muller"*  iFM),  is  already  quite  an  accomplish-  ilA2Rg  i,  with  7?^  being  the  universal  gas  constant.'**  For 

ment.  The  few  numbers  taken  from  Van  Vechten’s  calcu-  an  alloy  to  be  miscible  throughout  the  whole 

lations’’  (VV)  indicate  that  the  dielectric  model  predicts  concentration  range,  the  requirement  is  that  both  the 

results  at  larger  variance  with  experiments.  There  are  two  melting  temperatures  T,  and  F;  of  the  pure  4C  and  BC 

'  important  implications  about  the  FPT  that  can  be  drawn  compounds  be  greater  than  F,.  Table  VI  lists  the  values 

I  from  Table  III.  First,  the  theory  predicts  a  small  but  neg-  of  F^  associated  with  the  fl  values  in  the  FPT,  the  ratios 

ative  II  value  for  several  alloys.  This  not  only  means  that  7^ /F,  and  Tc/7';.  and  the  average  absolute  values  of  b;, 

there  is  no  miscibility  gap  in  these  alloys  but  also  implies  for  the  alloys  considered  in  Table  V.  In  Table  VT,  J ^  is 

!  a  tendency  toward  ordering,  in  which  the  substitutional  set  equal  to  zero  if  fl  is  negative  and  /  ,  is  chosen  to  be 

atoms  tend  to  be  surrounded  by  different  second-neighbor  the  lower  value  of  the  two  melting  temperatures,  so  the 

species.  For  stoichtometnc  compositions,  this  implies  a  criterion  for  not  having  .i  miscibility  cap  is  F,  <  I 
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There  is  an  empirical  rule^*’  stating  that  this  will  be  satis¬ 
fied  if  the  lattice  mismatch  j  Sq  |  between  the  two  alloy 
components  is  less  than  1.5%.  However,  we  find  that  (see 
Appendix  C)  a  more  precise  rule  is  |6(,  |  where 

6„  =  1.63V„  and  is  the  ratio  of  the  rms  bond-length 
amplitude  fluctuation  to  the  bond  length  at  the  melting 
temperature  T^.  The  values  of  T„  for  the  compounds  in¬ 
volved”  and  the  associate  values  estimated  from  Eqs. 
IC2)  and  IC3)  are  tabulated  in  Table  1.  The  model  used  in 
Appendix  C  yields  Tc/T^  =  (60^6^  *'•  This  suggests  that 
it  is  instructive  to  plot  Tc/T^  as  a  function  of  ,  60  '  /6„, 
as  is  shown  in  Fig.  4  for  the  calculated  from  FPT. 
This  plot  is  similar  to  the  A£  versus  60  curve  in  Fig.  2  be¬ 
cause,  in  fact,  fl  is  proportional  to  the  sum  of  the  A£ 
values  of  the  two  constituents  (see  Eq.  (41)].  However,  if 
Tc/T^  IS  plotted  against  i  alone,  the  FPT  points  are 
much  more  scattered,  and  those  of  SSHS  would  not  even 
exhibit  a  smooth  simple  quadratic  form  because  the  lower 
melting  temperature  Tj  is  not  a  smooth  function  of  8,1  ■ 
This  result  suggests  that  |  5o  /h„  <  1  is  a  better  cntenon 
than  I  80  I  <0.075.  Figure  3  also  clearly  shows  the  chem¬ 
ical  effects:  all  the  cation-suhstitution  alloy  points  lie 
below  the  solid  curve  and  all  the  anion-suhstitution  alloys 
have  [Tf/Ti)  values  on  or  above  the  curve,  corresponding 
to  negative  and  positive  shifts  in  IE  due  to  the  chemical 
terms.  Again,  the  curve  based  on  the  SSHS  model  is  an 
excellent  universal  representation.  From  the  figure,  we 
see  that  all  iP,Sb)  alloys  should  have  miscibility  gaps  and 
all  (As.Sb'  alloys  are  predicted  to  be  miscible,  although  on 
the  borderline,  because  the  actual  mixing  enihalpies  an 
larger  than  these  ideal  calculated  values  The  figure  also 
shows  that  ZniS.Te)  has  a  miscibility  gap  but  a  smalk-t 
value  of  Tf/T":  than  the  P.Sb  alloys,  despite  the  fact 
that  its  8(1 1  value  is  larger  All  these  previictions  are 
consistent  v, ith  the  available  experimental  evidence 

Vll,  SUMMARY  AND  CONCT-USUyN 

In  this  paper  a  simple  thev'rv  of  defect  substitulion  en 
ergies  is  formulated.  The  suhsutution  energy  is  compacih. 
separated  into  a  replacement  energy  and  a  distc'riion 
energy  of  the  pure  host  crystal  [see  Tq  2  1  However,  .1 
rigorous  application  of  this  theory  requires  an  improve¬ 
ment  in  certain  quantitative  tispects  of  Harrison's  bvinding 
theory,  p.irticularlv  the  elastic  constants  '  The  most  in¬ 
teresting  application  of  this  theory  presented  in  this  paper 
IS  its  perturbation  form  which  enables  us  to  absorb  the 
measured  elastic  coefficients  into  the  calculation  and. 
more  importantly,  to  study  the  chemical  effects.  The  ori¬ 
gin  of  chemical  infiuences  on  impurity  bond  relaxation 
can  be  attributed  to  three  mechanisms  [see  Eq.  il7)|  a 
chemical  force  £,,,  that  either  helps  or  hinders  lattice  re¬ 
laxation.  depending  on  whether  it  has  the  same  or  opp<i- 
site  sign  from  the  N-nui-length  difference  d-di  between 
the  host  and  impurity,  a  chemical  energy  that  depends  on 
the  diflerence  of  the  polanties  between  the  impurity  and 
host  bonds.  AC,,  and  an  effective  clastic  force  constant  H 
that,  when  positive,  also  tends  to  restrain  the  lattice  from 
ihstortioii  To  study  the  effect  of  boundanes  between  the 
core  atoms  around  an  impunlv  and  the  rest  of  the  clastic 
medium,  various  models  based  on  the  valence  force  field"' 


are  derived  and  their  results  are  compared  with  the  full 
perturbation  theory  and  available  experimental  data.  We 
found  at  least  five  models,  including  the  FPT,  that  pro¬ 
duce  the  correct  impurity  bond  lengths  with  variances  for 
the  compounds  studied  about  equal  to  the  experimental 
uncertainties  in  EXAFS'  *  '*  (~0.0l  .Al.  However,  some 
models  are  oversimplified  and  will  certainly  not  predict 
other  properties  equally  well.  However,  more  experimen¬ 
tal  lattice  constant  measurements  to  further  test  the 
theory,  particularly  on  (As,Sbl  and  (P,Sb)  substitution  sys¬ 
tems  for  which  there  are  larger  differences  between  dif¬ 
ferent  models,  are  needed  It  would  also  be  instructive  to 
sec  if  the  predicted  reversal  for  (Hg.CdiTe  is  found. 

The  excess  energies  of  impurity  substitution  are  also 
shown  to  prov  ide  good  c*slimates  of  the  mixing  enthalpies 
fl  of  pseudobinary  alloys.  The  chemical  shifts  are  found 
to  have  a  negative  net  coninbution  to  O  for  most  cation 
substitutions,  but  [xtsiiue  contributions  for  anion  substi¬ 
tutions.  The  chemical  reduction  of  !l  in  Ga.Inl  alloys  is 
larger  than  I  seal,  mole  ..^0  Ult.i'';  Several  VFF 
mexlels  and  the  full  perturbation  theory  prcduce  results 
for  n  that  are  as  gixxj  as  the  best  theories  with  one  adjust¬ 
able  paiamelei  Howevei.  the  full  thesuy  tends  to  yield 
answers  on  the  low  side  of  the  expenmenta!  values,  which 
we  argue  is  js  n  sfiould  be  bcsjuse  there  are  nonideal 
structures  ih.ii  a!s,'  contribute  to  !!,  I  he  calculated  fl 
values  and  tin  melting  ieniper.iiiiri.-s  aie  used  to  predict 
the  cxisternc  t  all  -v  inissibilitv  g.ips.  and  the  results 
correlate  wel.  wnh  ex;  eriments 

foiallv,  we  vs;^b  to  ^  .mmeiil  -’n  the  .K,ui.icy  of  the 
the-  lie'  tti.ii  .I't  .  '.lie,  led  li’  the  pieseiil  in- del,  I  he  per- 
lurb.tii.-n  the  '  .'i.i'  .lite.Klv  t-H-eii  streuhed  beyond  its  ex- 

pe..-evi  legi  -n  -t  '.i.i.ti'v  jiid  piedkls  t-'  within  expen- 
tnent.i.'  .oki  rt.iiMiev  -  o,  .y  even  t.-r  ,.tses  with  large 
b-'ii  l  .engit;  .ntterei.  t-s  -  '  I  I  hi  i, ,  ur.icy  can  only 
N-  itiipi  veil  .t  tie.  tui!  n.-Tiper-urbation  theory  outlined  in 
Sfs  II  .'  .iM-,l  Ih.s  .. .:i' uiat U'-i  IS  neeOvst  for  the  strong 
subsiniiti-'ti  in.it  w  re  r,-n  ,  'iisutered  in  this  paper 

exainpies  ii,  li  ( la  H  In  ,  N,P  .  N.Sb  ,  O.S'. 

■(t.Se  .KHi  <)  le  -utv'i  it  111  loii'  Aliln-ucn  we  tseiieve  that 
lor  the  proj  •criii-'.  I'e.ited.  tne  in, ski  with  a  c  nlinuuni  at- 
t. Is  heel  to  ihe  sc,,  tut  'hed  is  .is  .i.'urati  as  the  perturba¬ 
tion  ihes'iy  use'tl.  11  reni.iins  te'  be  se-cn  if  this  is  true  fcir 
other  pr.'perties.  espe,  lally  strain  '-'elTkients  hinally,  the 
present  thes'ry  h,is  been  extended  t,-  study  alloys'*'  by 
enibediiing  clusters  in  an  effivtive  medium  This  enables 
us  tsi  stuily  the  bond  length  and  energy  sanations 
throughout  the  whole  concentration  range  However,  a 
quantitative  calculation  still  aw.iits  an  improvement  of  the 
accuracy  of  Harrison's  theory  .-A  similar  priKedure  is 
also  being  extendevi  to  a  study  of  the  alloy  electronic 
structure  for  which  a  cluster  Cl’.-\  coherent -potential  ap¬ 
proximation'  involving  both  pt'tential  and  structural  disor¬ 
der  '  will  be  used 
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APPENDIX  A:  ELASTIC  ENERGY  IN  CONTINUUM 


In  Sec.  Ill  the  elastic  energy  outside  a  sphere  of  radius 
R  centertxl  at  the  impurity  is  assumed  to  be  a  continuum 
with  a  radial  displacement  uac?/r^.  If  the  displacement 
at  R  is  ug,  then  u{i)  =  Uo{R'/r‘'lT.  The  energy  density  in 
the  continuum  is  given  by 

betrl  =  |C|i<eij  -rCn  I 

4"  ■;  C44  (  -t-  6y2  "h  1  , 


2.  a  term.s  from  the  .second-shell  bonds,  ^  terms 
between  the  first-  and  second-shell  bonds 
and  among  the  second-shell  bonds 

For  these  terms  we  need  to  consider  the  four  bond  vec¬ 
tors  pointing  away  from  C.  They  are 

r|=(  -1  4-8.  -1+8,  -l+8W/v'3  . 

r2  =  (-  l+8,  1+8  +  y,  1  +  6  +  y)d/cl  , 

r3  =  (l+8  +  y,  — 1+8,  1  +y  +  8W/v^, .... 

Then  A(r2-r2)=  ft8  +  2y A(r|-r,)=  -  |(8  +  y)£/^,  and 
A(r2  r))  =  48d^.  Thus,  the  a  terms  from  the  second-shell 
bonds  become 

4x3x3a[T(8  +  2yW^]V8d^  =  2a(8  +  28)^d’  . 


where 


bu. 


e„  =  ----  =R^Uoir'  —  3x^)/r-  , 

dx 


duj,  duy 

3,  3x 


=  -— +  — =  -6R^Uoxy/r 


the  P  term  between  the  first-  and  second-shell  bonds  are 
4x3x3)9(y<8+rW^]V8d’^  =  2(r+3)‘iftf^  , 
and  the  p  terms  among  the  second-shell  bonds  are 
4x3x3^(T8d^)V8d^  =  28^j3d^  . 


Thus,  the  total  elastic  energy  in  the  continuum  is 

=  4nRug(  tC||  —  yCi2  +  ■JC44) 

I  =CRul  , 

where  the  effective  shear  coefficient  is  given  by 
.  C  =  Trt  1 .  bC'  11—1.  bC*  1 2  ^  *4  ^  • 

I 


3.  a  terms  for  the  third-shell  bonds,  terms 
between  the  second-  and  third-shell  bonds 
and  among  the  third-shell  bonds  adjacent 
to  the  second-shell  atoms 

For  these  terms  we  need  to  consider  the  bond  vectors 
pointing  away  from  B  in  Fig.  1.  They  are 

r2  =  (  1  -6,  -  1  -  y  —  6,  —  1  —  y  -  bid /v'l  , 

ri  =  (l  +  y',  l  +  3y'-y,  l  +  3y'-y)d/v^  , 

rj  =  (  -  1  -y”,  1  +3y''-y,  —  1  +  y”  -y)d /t/3  , 


APPENDIX  B:  DISTORTION  ENERGY 

In  this  Appendix  we  count  the  detailed  contributions  of 
the  bond-stretching  terms  A(r,  r,  I  and  “bond-bending" 
terms  Atri'r^l  for  i-r  j  in  VFF  [Eq.  tl3))  that  enter  Eq. 
(16)  in  FPT  and  in  the  VFF  models  in  Sec.  III. 

I.  a  and  B  terms  from  the  first-shell  bonds 

The  four  bond  vectors  pointing  away  from  the  central 
impurity  according  to  Fig.  1  are 

T)  -  (  1  b,  1  b,  1  b  id  /  \  . 

;  r-'^i  1  b,  -  1  -  b.  1  bidv  . 

Thus,  AGi  Ti)  2.Vi‘  .inJ  A'r,  r;  The  u 

terms  contribute 

I  4'-  3n/‘  -2bd‘i'/8d'  iS<(;b'd'  , 

and  the  P  terms  contribute 

b  •  3/.</'  ‘  8d  'Hd-  /3;b'd-  . 

If  an  A  atom  is  replaced  b\  a  /?  atom,  as  was  tione  in 
FPT.  the  (t I  and  Pi  are  replaceil  bv  ti  ami  /)’.  rcspcvtivelv 

I 


and 

r,  =  (  - 1  -  y",  —  1  -I-  y "  —  y ,  1  i-  3y "  -  y  Id  V  3 
Thus,  we  have 


A(r,T2)  = 

-  ^136  + 

Sr') . 

Airor  ,!  -  - 

Yt6-2y 

-5y”) 

At  rir4) 

jd^-"  . 

A:  f  ;  r,  1 

1  3, 

At  r,r,i 

tbs” I  , 

Air|-r|i 

4  Ml,- 

4y)  . 

model  .A, 

-  r".  0, 

so  the  I 

and 


Kinds  K’conie 

4  ■  3  ■  3(([  2  Ar .  r,  •  '  Ar,  +1  ’  j  8d’  S.iiy'd’  . 
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the  0  terms  between  the  second-  and  third-shell  bonds  are 
4x3x3/3l(Ar|T2)^  +  2(Ar,-r3)^]/8d^ 

=  \{b  +  2Y)^]=  TPl>^d^  +  0(8  +  2Y)^d^  , 


f)  =  (  -  1  -r  y,  -  I  -  y,  -  1  )d/v'i  , 

r2=(-l,  I,  lW/t/3  . 

rj  =  (  1,  -  1  -i-y,  1  -t-y)d/vl  , 


and  the  0  terms  among  the  third-shell  bonds  adjacent  to 
the  second-shell  atoms  are 


4x3x3/3(2(Ar,T3)V(8d*)-|-(Ar3T4)^]=40v’d-  . 


For  continuum  and  the  only  contribution  from  this  group 
are  the  0  terms  between  the  second-  and  third-shell  bonds. 
Since  the  displacements  in  the  continuum  are_proporti^al 
to  l//^^y'  =  8v^r/(19v'l9)  and  y"  =  8t/2y/( llv/ll). 
Thus,  these  0  terms  become 


y  ( 36  +  5y' ^ (6 2y  -  5y'' 


=  7/3 


40t/2 

36 -f - 

19v/T9 


-1-2 


C  ,  ^  40v^ 

TI7Tr>- 


rf'. 


4.  0  terms  for  bonds  adjacent 
to  the  third-shell  atoms 

These  terms  only  enter  model  A,  so  r'  =  r"  =  0.  There 
are  two  different  groups,  one  like  those  adjacent  to  C'  and 
another  like  those  meeting  at  C".  The  four  bond  vectors 
pointing  away  from  C'  are 

fi  =  (  -  1,  -  1  -f  y,  -  1  -(-y  W/v  3  , 

r2  =  (-l,  1.  IW/v'l  , 

r,  =  (l,  -1,  IW/v'l  , 

and 

r4  =  (I,  1.  -l)d/v'3  . 

Thus,  the  only  contribution  from  this  group  is 
4x3x30tA(r-r2)]^/8d^  =  2/3y'd’  . 

The  four  bond  vectors  around  C"  are 


and 

r4  =  (l,  1,  -l)d/v^  , 

which  only  results  in  the  first-order  term  A(r2T3)=  yyd^. 
Thus,  the  group  contributes  to 

4v3x3^[2(Ar2  r,)^]/8d^  =  40y^d*  , 

and  the  combined  contribution  from  these  two  groups  is 
b0Y^d\ 

APPENDIX  C:  CRITERION  OF 
MISCIBILITY 

Starting  with  Eq.  (24)  and  using  the  SSHS  model 
a— at,  one  finds  the  mixing  enthalpy  parameter  12  to  be 

n  =  ya(d4c-rf8r>%  .  (Cl) 

where  .Vq  is  Avagodro’s  number  and 

Then  relate  the  mean-square  bond-length  fluctuation 

at  the  melting  temperature  T„  to  T„  for  a  compound  by 

equating  the  average  potential  energy  per  unit  cell  to  half 

of  the  thermal  energy: 

<  I%,>-.4XTa(f-)=4(2x3fc8r„)  ,  (C2) 

where  kg  is  the  Boltzmann  constant.  Defining  a  Lieder- 
mann  ratio  of  melting  X„  by 

(<;’))’"’;=  V„d  (C3) 

and  choosing  the  mixing  criterion  to  be  T^/T„  <  1,  where 
T„  now  is  the  smaller  value  of  the  two  melting  tempera¬ 
tures  of  the  constituent  compounds,  we  require  that 

34  flk  3  (t^.4C 

_ - ^ ^  ^ — - ^  1  ((^4) 

T„  AR^aX'„d^  8  X-„d- 

or 

I  60  V6„  <  1  ,  (C5) 

where  6„  =  1.63V„,  and  60  is  the  percentage  bond- 
length  difference. 
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The  sensitivity  of  defect  energy  levels  in  semiconductors  to  the  host  band  structures  and  impurity 
potentials  has  been  studied  for  approximately  30  impurities  in  CdTe  using  four  different  band- 
structure  models.  The  discrepancies  in  the  defect  levels  between  two  different  sets  of  band  struc¬ 
tures  and  impurity  potentials  are  found  to  range  from  less  than  0.1  eV  to  the  whole  band  gap  (1.6 
eV).  The  band-structure  effects  are  analy2ed  here  in  terms  of  detailed  partial  densities  of  states. 
Examples  of  contradictory  predictions  from  different  band  structures  are  illustrated,  and  ways  to 
improve  the  theory  are  suggested. 


I.  INTRODUCTION 

In  several  of  our  recent  papers,’"^  we  have  applied  a 
method  to  calculate  the  band  structure  of  semiconductors 
that  is  both  efficient  and  accurate.  Because  the  procedure 
involves  casting  the  basis  functions  into  orthonotmal  local 
orbitals*  (OLO),  our  method  has  the  advantages  common 
to  empirical  tight-binding  (ETB)  calculations,*”'®  except 
that  the  Hamiltonian  matrix  elements  to  all  ranges  are  re¬ 
tained.  The  inclusion  of  these  higher  coefficients  makes  it 
possible  to  produce  excellent  band  structures  including 
conduction  bands  and  effective  masses.  The  method  also 
yields  wave  functions  for  optical  property  calculations.’ 
Moreover,  its  OLO  description  also  permits  its  extension, 
through  the  coherent-potential  approximation,  to  al¬ 
loys,’”’ 

The  recent  attention  focused  on  defects  in  semiconduc¬ 
tors  motivated  us  to  apply  our  method  to  this  problem. 
The  theories  of  defects  have  ranged  from  very  sophisticat¬ 
ed  self-consistent  density-functional  theory"””  (SCDF) 
to  simple  ETB  calculations.  It  is  generally  recognized 
that  SCDF  is  as  accurate  in  defects  for  the  ground-state 
properties  as  it  is  for  pure  semiconductors,  but  less  certain 
in  assigning  excited  energy  levels.  ETB,  because  it  can 
produce  results  for  many  systems  in  one  study,  claims  to 
predict  the  trends  of  deep  levels'®  even  if  the  accuracy  for 
a  given  impurity  may  be  poor.  However,  this  contention 
remains  to  be  verified. 

To  assess  this  concern,  we  ask  the  following  question; 
“How  sensitive  are  defect  levels  to  host  band  structures 
and  impurity  potentials?"  To  this  end,  we  have  adopted 
the  simple  yet  nontrivial  defect  model,  that  of  site- 
diagonal  substitutional  defects  often  used  in  ETB  studies. 
CdTe  was  selected  in  this  study  because  its  band  structure 
has  been  examined  in  great  detail  by  us,  and  there  are 
three  published  band-structure  models*"'”  that  we  could 
easily  generate  for  comparison.  There  is  also  a  consider¬ 
able  body  of  experimental  data  on  deep  states  in  this  sys¬ 
tem.'*”'’ 

.11 


U.  CALCULATIONAL  PROCEDURE 

In  the  simple  site-diagonal  substitutional  defect  model, 
the  impurity  energy  levels  E  are  determined  by  the  equa¬ 
tion 

l-u„g„(£)=0,  (1) 

where  a  designates  the  symmetry  of  a  local  state,  e.g.,  P j, 
r 7,  and  Pg  on  an  atomic  site  in  the  zinc-blende  structure, 
and  go  is  the  real  part  of  the  diagonal  matrix  element  of 
the  host-crystal  Green  function,  go  can  be  calculated 
from  the  partial  density  of  states  (PDOS)  by 

ggiE)-  J pa{e)/{E  ~f)de  .  (2) 

The  PDOS  is  given  by 

Po(f)  =  2  !«“<*)  I  .  (3) 

n.k 

where  k)  are  band  energies  and  a°{k)  are  the  probabili¬ 
ty  amplitudes  of  the  band  state  in  the  Bloch  basis  con¬ 
structed  from  the  OLO  labeled  by  a.  The  Brillouin-zone 
integration  in  Eq.  (3)  is  calculated  using  an  accurate  ray 
scheme.'* 

Because  a  principal  concern  of  this  paper  is  the  sensi¬ 
tivity  of  impurity  levels  to  the  host  band  structures,  we 
should  emphasize  the  difference  between  our  method  and 
ETB.  Our  method  consists  of  four  steps. 

(1)  We  start  with  four  Gaussian  orbitals  per  atom  and 
empirical  pseudopotentials,'’  and  compute  the  Hamiltoni¬ 
an  matrix  Hik)  and  overlap  matrix  Sik)  as  was  done  by 
Kane’®  and  Chadi.’' 

(2)  The  Gaussian  orbitals  are  transformed  into  OLO,* 
so  //(/()  is  transformed  into  Ho^k)  and  5  into  the  identity 
matrix.  The  band  structures  calculated  from  Hglk)  are 
accurate  to  5%  as  compared  to  more  sophisticated 
methods  using  the  same  potential.' 

(3)  A  spin-orbit  Hamiltonian  in  the  OLO  basis*  is  in¬ 
corporated  to  deal  with  this  interaction. 

(4)  To  compensate  for  the  effects  of  truncated  basis  and 
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nonlocal  potentials,  a  perturbation  Hamiltonian  is 
added.  H  \  has  the  same  form  as  a  truncated  ETB  Hamil¬ 
tonian.  The  parameters  in  H\  are  adjusted  to  fine  tunc 
the  important  band  energies  and  effective  masses. 

Although  both  ETB  and  our  methods  are  empirical, 
there  are  two  major  differences. 

(1)  While  most  ETB  retains  the  H  matrix  elements  only 
to  the  first-  or  second-neighbor  shell,  ours  extends  to  all 
ranges,  so  that  the  high  Fourier  components  needed  to 
produce  the  sharp  band  curvatures  are  properly  given. 

(2)  Our  method  can  directly  generate  wave  functions  for 
calculation  of  other  properties. 

Thus,  while  our  method  yields  more  accurate  band  struc¬ 
tures,  it  retains  much  of  the  advantage  of  ETB,  namely 
the  computational  speed  and  a  simple  direct-space 
description  of  the  Hamiltonian. 

111.  BAND  STRUCTURES  AND  PARTIAL 
DENSITIES  OF  STATES 

Figure  1  depicts  the  four  band  structures  to  be  con¬ 
sidered  for  CdTe.  Our  result  is  in  panel  (a);  panels  (b) 
(Ref.  8)  and  (c)  (Ref.  9)  are  two  ETB  band  structures  with 
the  Hamiltonian  matrix  elements  truncated  at  second 
neighbors.  (Because  different  parameters  were  selected, 
these  two  band  structures  are  not  identical.)  Panel  <d) 
(Ref.  10)  results  from  the  use  of  five  basis  orbitals  per 
atom;  the  extra  one  is  an  excited  s  state.  All  these  band 
structures  are  adjusted  to  have  the  proper  fundamental 
band  gap  of  1.6  eV.  The  principal  differences  one  sees  on 
first  inspection  are  in  the  band  curvatures,  especially  the 
conduction  bands.  The  effective  mass  at  the  bottom  of 
the  conduction  band  in  panel  (a)  is  0.1  times  the  free- 
electron  mass,  in  agreement  with  experiment,”  while  in 
other  panels  it  is  more  than  twice  as  large. 

Figure  2  shows  the  densities  of  states  (DOS)  for  each  of 
the  band  structures  in  Fig.  1.  While  the  valence  bands  at 
least  exhibit  general  common  features,  the  conduction 
bands  are  almost  unrecognizable  as  representing  the  same 
compound.  In  panels  (c)  and  (d),  for  example,  there  is  a 
second  band  gap  above  the  fundamental  gap.  Also  note 
that  there  are  two  extra  narrow  peaks  associated  with  the 
two  extra  excited  s  orbitals  (one  for  Cd  and  the  other  for 
Te)  included  in  the  calculation. 

To  analyze  the  band  effects  on  defect  levels  [see  Eqs.  (!' 
and  (2)],  the  DOS  is  further  decomposed  into  partial  den¬ 
sities  of  states  for  F^ts),  Filp'^’),  and  states  on 

the  Cd  and  Te  sites,  as  shown  in  Figs.  3—6.  The  F^ 
PIX)S  are  not  shown  because  they  are  nearly  the  same  as 
Ft  with  only  a  slight  upward  energy  shift.  These  PIX)S 
show  how  the  ‘'atomic"  levels  evolve  into  hand  states. 
These  curves  contain  useful  information  about  many 
properties,  e  g.,  the  relation  between  the  crystal  bonding 
and  atomic  energies,  and  how  potential  disorder  in  alloys 
affects  different  pans  of  the  bands,’ ' '  in  addition  to  de¬ 
fect  levels  studied  here. 

The  F^tCdi  PIX3S  shown  in  Fig.  3  split  between  the 
conduction  and  valence  bands.  It  is  generally  assumed 
that  the  cation  s  stales  in  lll-V  and  ll-Vl  compounds 


FIG.  1.  Four  band  structures  of  CdTe  used  for  comparative 
studies:  (a)  present  work,  (b)  Ref.  8,  (c)  Ref.  9,  and  (d)  Ref.  10. 
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FIG  ^  Densities  of  slates  calculated  from  the  four  band 
structures  in  Fig  1 
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evolve  into  the  conduction  bands,  while  the  anion  p  states 
make  up  most  of  the  major  valence  bands  just  below  the 
gap.  Thus  it  is  perhaps  a  surprise  to  see  a  prominent  peak 
denved  from  the  cation  s  states  at  the  bottom  of  the  ma¬ 
jor  valence-band  structure.  However,  this  is  a  general 
feature  for  all  sp ’-based  compound  semiconductors. 
These  are  the  states  responsible  for  the  first  observed 
breakdown  of  the  virtual-crystal  approximation  for  a 
semiconductor  alloy;  Hg|_;,Cd,Te  (which  is  caused  by 
the  large  j-energy  shift  between  the  Cd  and  Hg  sites).*'*'^^ 

A  more  detailed  examination  draws  attention  to  some 
important  differences  among  the  four  panels  in  Fig.  3; 
the  valence-band  peak  in  panel  (c)  is  about  2  eV  higher 
than  the  rest,  and  it  is  also  high  compared  to  experi¬ 
ment.’^  Our  conduct  ion -band  PDOS  in  panel  (a)  is 
broader  than  the  others.  The  ratio  of  the  integrated 
PDOS  in  the  conduction  bands  to  that  in  the  valence 
bands  in  our  model  is  larger  than  those  in  other  panels. 
Also  our  PDOS  just  below  the  valence-band  edge  is  obvi¬ 
ously  smaller  than  that  found  in  other  models. 

Figure  4  shows  that  the  Cd  p  states  are  concentrated  in 
the  conduction-band  states.  This  is  particularly  true  in 
panel  (a),  where  their  contribution  to  the  valence-band 
states  shrinks  almost  to  nothing.  In  other  panels,  there 
are  still  sizable  (  —  20%)  valence-band  states.  In  contrast, 


all  four  panels  in  Fig.  5  show  that  the  Te  s  states  are  con¬ 
fined  to  the  deep  valence-band  states,  as  generally  recog¬ 
nized.  Finally,  Fig.  6  shows  that  the  Te  p  states  dominate 
the  upper  valence-band  states.  Panel  (a)  has  much  less 
conduction-band  content  than  the  other  three  panels.  As 
we  will  see,  these  differences  can  result  in  quantitatively 
or  even  qualitatively  different  predictions  about  the  deep 
levels. 

IV  IMPURITY-LEVEL  DETERMINATION 

A  convenient  way  to  study  the  impurity  energy  levels 
using  Eq.  (1)  is  to  rewrite  it  as  Cq  =  1  /g^i  E)  and  plot  E  as 
a  function  of  v.  Once  this  E-v  curve  is  deduced  for  each 
a.  the  deep  levels  E^  for  a  given  impurity  can  be  read  off 
the  curve  by  drawing  a  vertical  line  at  the  appropriate 
value  of  i’„  for  the  impurity.  We  set  the  zero  of  energy  at 
the  top  of  the  valence  bands.  Because  the  gap  is  1.6  eV, 
we  will  focus  on  levels  in  the  energy  range  from  —0.5  to 
2.0  eV. 

Calculations  have  been  performed  for  all  neutral  impur¬ 
ities  listed  in  Table  I.  Because  we  do  not  believe  that 
there  exists  a  uniformly  accepted  table  for  c  we  have 
adopted  a  table  that  we  used  for  structural  studies.”’’* 
Table  1  lists  the  term  values,  which  we  obtained  from  to- 


TABLE  I.  j-  and  p-state  correlated  term  values  in  units  of  -eV  The  top  entry  is  the  i-state,  the 
second  the  pi/2-state,  and  the  third  the  pj/j-stale  energy.  (All  energies  are  negalive.) 
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FIG.  7.  The  E-u  curves  for  the  Ft  states  on  a  Cd  site. 

tal  energy  differences  between  atomic  configurations  cal¬ 
culated  using  the  norm-conserved  pseudopotentials^’  and 
self-consistent  charge-density-functional  theory,  with  the 
first  ionization  energies  adjusted  to  be  the  experimental 
values.^*  These  term  values  are  found  to  yield  consistent¬ 
ly  better  structural  properties^’  in  Harrison’s  theory’’’^* 
than  those  based  on  Mann's  values’^  adopted  by  Har¬ 
rison.^*  The  impurity-potential  parameters  will  then  be 
taken  as  the  difference  of  the  term  values  between  the  im¬ 
purity  atom  and  Cd  <or  Tel.  To  study  the  sensitivity  of 
Ea  to  Va,  we  shift  I’a  by  ±0.5  eV  and  compute  the  corre¬ 
sponding  changes  in  the  energy  levels. 

Figures  7—10  display  the  E-v  curves  for  several  a. 
Each  figure  has  four  curves,  corresponding  to  the  four 
panels  of  PDOS  in  each  of  Figs.  3-6.  The  functional 
behavior  of  these  curves  can  be  understood  qualitatively 
using  Eq.  (2)  and  Figs.  3—6.  If  E  lies  in  the  gap,  the  con¬ 
tribution  from  conduction  bands  is  negative,  but  positive 
from  the  valence  bands.  The  closer  the  PIX)S  to  the  E  in 
question,  the  larger  will  be  its  influence.  Applying  this 
argument  to  the  FjiCd)  representation,  we  see  that  the 
curves  in  Fig.  7  are  negative  in  the  gap  region  because  the 
PDOS  in  Fig.  3  near  the  bottom  of  the  conduction  bands 
are  much  larger  than  those  near  the  valence-band  top. 
Thus,  on  the  Cd  site,  only  impurities  with  an  j  energy 


C  Te 


-25  -20  -15  -10  -5  0  5  10  1  5  20  25 

V  (eVI 

f-TG  X.  The  £-1'  curves  for  ihc  T -  states  on  a  Cd  site 


N*  Rb 
Cd  Cu  Li  I  I  Ci 


5  to  15  20  25  30 

V  lev  I 


FIG.  9  The  E-v  curves  for  the  F^  states  on  a  Te  site. 

below  the  Cd  s  level  ( —  8.99  eV)  will  produce  a  r6  level  in 
the  gap.  However,  we  note  that  in  Fig.  7,  g„(£)=0  for 
models  (b)  and  (c)  just  below  the  valence-band  edge  be¬ 
cause  of  cancellation  between  the  conduction-  and 
valence-band  contributions.  At  this  E  value,  the  E-v 
curve  switches  from  v—  —  oa  to  v=ao  (not  shown);  an 
ideal  vacancy  level  (corresponding  to  Cq  =  oo )  is  located  at 
this  E.  A  similar  consideration,  but  with  the  conduction 
and  valence  bands  interchanged,  leads  to  an  understanding 
of  the  curves  in  Fig.  10.  Using  the  same  principle,  we  can 
easily  understand  why  all  curves  in  Fig.  9  for  the  TjITe) 
representation  are  positive,  but  the  reasons  for  the  large 
displacements  between  these  curves  are  not  easy  to 
deduce.  In  Fig.  8,  the  curve  labeled  a  is  distinctively  dif¬ 
ferent  from  other  curves,  because  the  PDOS  in  panel  (a)  in 
Fig.  4  is  completely  dominated  by  the  conduction  band; 
however,  for  the  other  panels  the  PDOS  just  below  the 
valence-band  edge  are  as  large  as  those  just  above  the 
conduction-band  edge.  This  produces  a  very  sharp  nega¬ 
tive  E-v  curve  for  (a),  but  split  behavior  for  (b),  (c),  and 
(d). 

These  E-v  curves  provide  a  clear  picture  of  how  dif¬ 
ferent  host  band  structures  may  affect  the  deep  levels. 
Numerical  values  for  the  impurity  levels  can  be  obtained 
from  these  figures  by  drawing  vertical  lines  at  the  ap- 
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propriate  impurity  potentials  (i.e.,  differences  between  the 
term  values  listed  in  Table  I),  as  has  been  shown  for 
several  represertlative  impunlies.  To  provide  a  more 
quantitative  comparison.  Table  II  lists  some  calculated 
impurity  levels  £„  and  the  corresponding  changes  ikE^ 
due  to  the  l-eV  change  in  v^- 

V  RESULTS  AND  CONCLUSION 

To  summarize  we  recall  that  band  models  (b)  and  (cl 
are  the  same  second-neighbor  ETB  with  two  different  sets 
of  parameters,  and  model  (d)  is  a  first-neighbor  ETB  with 
one  extra  s  orbital  per  atom.  Our  model  (model  (a)]  has 
the  form  of  ETB  but  is  derived  in  a  very  different  manner 
and  includes  all  the  long-range  interactions.  Therefore, 
we  expect  that  the  results  from  models  (b)  and  (c)  will  be 
close,  model  (d)  will  have  larger  discrepancies  from  (b) 
and  (c)  than  that  between  (b)  and  (c),  and  model  (a)  will 
differ  even  more.  This  is  evident  from  Figs.  7—10  and 
Table  II.  We  found  the  energies  for  the  r,(Cd),  r7(Te), 
and  FglTe)  states  produced  by  models  (bl  and  (c)  agree 
within  0.1  eV.  For  the  other  states,  i.e.,  FtlTe),  r7(Cd), 
and  r glCd),  the  energies  from  (b)  and  (c)  are  qualitatively 
similar,  but  the  difference  can  be  as  large  as  0.4  eV.  The 
largest  discrepancy  between  models  (d)  and  (b)  [or  (c)]  is 
more  than  0.5  eV,  and  that  between  (a)  and  other  models 
is  more  than  1  eV.  The  largest  difference  comes  from  the 


p  levels  on  a  Cd  site.  For  example,  the  filled  p  level  of  C 
on  a  Cd  site  in  model  (a)  is  a  resonance  state  just  below 
the  valence-band  edge  but  is  a  donor  state  in  the  other 
models.  Similarly,  model  (a)  puts  the  neutral  Te  antisite 
defect  p  levels  at  about  -j-  and  of  the  gap  [£(r7 1  =  0.48 
eV  and  £(rg)=0.95  eV],  while  other  models  assign  them 
as  resonance  states  inside  the  conduction  bands.  We  also 
note  that  the  discrepancies  between  different  models  are 
not  uniform,  but  vary  with  Consider  FftlCd)  for  ex¬ 
ample.  All  four  models  yield  the  same  ordering  and  about 
the  same  energies  for  the  group-III  impurities  Al,  In,  and 
Ga.  However,  as  i>  becomes  more  negative,  the  splitting 
between  the  curves  increases,  so  the  discrepancies  become 
larger  [  —  1  eV  difference  between  models  (a)  and  (d)  for  I 
impurity].  Similarly,  for  the  r7(Te)  states,  all  four  models 
put  the  Sn  impurity  energies  close  to  the  valence-band 
edge,  but  the  agreement  deteriorates  as  increases. 

Regarding  the  sensitivity  of  energy  levels  to  impurity 
potentials.  Table  II  shows  that  a  1-eV  shift  in  produces 
a  change  in  ranging  from  less  than  0.1  to  0.6S  eV. 
Very  little  is  known  about  the  size  or  trends  in  errors  in¬ 
troduced  in  t'o  from  the  use  of  atomic  term  values.  How¬ 
ever,  we  know  that  the  discrepancy  of  Vg  between  two  dif¬ 
ferent  tables  of  atomic  term  values  can  be  larger  than  2 
eV.  This  discrepancy  translates  into  an  uncertainty  of  less 
than  0.1  to  more  than  1  eV  in  the  impurity  energy  levels. 


TABLE  II 
units  of  eV .  F 

Defect  energy  levels  £  and  changes  iF  due  to 
i)  stands  for  ideal  vacancy. 

a  1-eV  change 

in  the  impurity-potential  parameter. 

All  energies  are  in 
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0.39 

1.42 

0.24 

1.33 

0.23 

1.57 
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which  is  comparable  to  that  due  to  different  host  band 
structures. 

Putting  this  large  uncertainty  in  the  deep  levels  against 
a  band  gap  of  1.6  eV,  we  are  left  with  great  doubts  about 
the  predictability  of  this  oversimplified  theory.  Unfor¬ 
tunately,  the  experimental  means  available  for  identifying 
micrixlefects  in  semictmductors  are  still  very  limited,  and 
the  ah  initio  band  theory  is  still  not  capable  of  accurately 
predicting  the  energy  levels.  Thus,  there  is  a  great  temp¬ 
tation  to  use  simple  theories  like  the  one  earned  out  here 
to  help  with  the  identifications.  To  illustrate  this  point, 
consider  the  following  examples:  Table  II  shows  that  Li 
on  a  Te  site  has  an  v  level  of  0  14  eV  in  model  (a),  so  one 
may  be  tempted  to  relate  it  to  the  acceptor  state  identified 
experimentally.''*  However,  this  is  not  the  hydrogenic  ac¬ 
ceptor  stale  on  a  Cd  site,  as  one  might  anticipate.  One 
might  also  want  K'  assign  the  7  and  7  g^P  states  for  the 
Te  antisite  p  levels  on  the  Cd  site  found  from  model  (a)  as 
those  seen  in  experiments.''  Because  of  the  large  uncer¬ 
tainty  in  the  calculation,  the’se  results  should  be  regarded 
as  suspicious  surpn.ses  rather  than  thetiretical  confirma¬ 
tions. 

The  results  presented  here  should  not  discourage  con¬ 
tinued  research  on  the  ETB  approach,  but  improvement  is 
clearly  needed.  Work  ranging  from  universal*''^’''**  to 
specific^*' structural  studies  to  our  band  calculations 
and  alloy  studies'"^  indicates  that  the  ETB  type  of  theory 
is  practical  for  both  bonding  properties  and  electronic 
structures.  The  reason  that  ETB  works  well  for  some 
properties,  e.g.,  photoemission  spectra  and  bonding  prop¬ 
erties,  but  not  for  impurity  levels,  is  that  the  former  de¬ 
pend  only  on  the  gross  total  density  of  slates,  while  the 


latter  have  been  shown  to  be  sensitive  to  the  details  of  the 
partial  densities  of  states. 

To  establish  the  credibility  of  ETB  in  defect  studies, 
one  needs  to  look  at  the  problem  more  seriously.  The 
most  difficult  and  yet  important  task  is  to  develop  a  belter 
way  for  determining  the  Hamiltonian  matrix  elements. 
Haas  et  al*  and  Harnson’''^"  have  suggested  using  the 
atomic  term  values  as  the  diagonal  matnx  elements.  Our 
work'~^  has  suggested  using  a  universal  long-range  in¬ 
teraction  to  improve  the  accuracy  of  the  conduction 
bands.  Several  studies''^''’*  '^  have  also  pointed  out  scal¬ 
ing  rules  of  the  matnx  elements.  A  combination  of  these 
ideas  may  lead  to  an  acceptable  model.  Secondly,  both 
the  bonding  and  deep-level  states  of  impunties  should  be 
studied  at  the  same  time  in  order  to  provide  correlated  in¬ 
formation  for  defect  identification.  Finally,  more  realistic 
models  should  be  examined.  Besides  the  substitutional 
site-diagonal  defects,  one  should  consider  the  possibility 
of  interstitial,  paired,  and  even  more  complex  defects. 
One  also  needs  to  deal  with  long-range  impunty  poten¬ 
tials,  possible  charge  shifts,  and  lattice  distortions.  Pro¬ 
gress  in  all  these  areas  can  be  expected  if  the  calculation  is 
constantly  correlated  with  experiments  and  available 
ab  initio  theory. 
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A  calculation  is  presenied  to  explain  the  anomalous  experimental  behavior  of  the  Si  2p  core- 
exciton  binding  energy  and  linewidth  in  Sij,Gei-,  alloys  The  observed  minimum  in  the  linewidth 
near  x  0.15  can  be  explained  as  the  result  of  a  competition  between  intrinsic  broadening  due  to 
screening  and  extrinsic  alloy  broadening.  For  pure  Si,  the  binding  energy  is  estimated  to  be 
0.15  ±0.05  eV  and  the  width  is  shown  to  be  smaller  than  that  observed  at  x  =  0.15. 


PACS  numbers:  7l,55.Fr,  71.35. +z,  78. 70.Dm 

Until  1984,  the  Si  2p  core  exciton  was  believed  to 
have  an  anomalously  large  binding  energy.'*'®  Later, 
Newman  and  Dow"  proposed  a  radically  different  pic¬ 
ture  in  which  the  St  2p  core  exciton  is  in  fact  a  reso¬ 
nance  with  a  negative  binding  energy.  They  further 
predicted  that  the  exciton  binding  energy  remains  neg¬ 
ative  throughout  most  of  the  Sij,Gei_;t  alloy  composi¬ 
tion  range,  except  near  .v  —  0.20  where  it  becomes 
positive.  In  a  recent  experiment,'^  Bunker  et  al.  found 
an  anomalous  sharpening  of  the  exciton  spectra  near 
x  =  0  15,  the  data  were  interpreted  to  support  the 
Newman-Dow  point  of  view.  Yet  the  most  recent  ex¬ 
periment'^  still  suggests  a  positive  value  for  the  bind¬ 
ing  energy  in  silicon. 

In  ihis  Letter,  we  present  a  calculation  that  offers  a 
plausible  resolution  to  the  above  problem.  In  our 
theory,  the  calculated  Si  2p  core-exciton  binding  ener¬ 
gy  £*(  v)  and  the  linewidth  A(x)  in  Si;,Gei_,  alloys 
are  sensitively  dependent  on  the  parameter  E®.  A 
comparison  of  the  calculated  A(x)  with  the  experi¬ 
ment"  suggests  a  positive  value  0.15  ±0.05  for  Ej. 
The  anomalous  experimental  spectrum'^  near  .x  =  0  1 5 
IS  explained  as  a  result  of  a  competition  between  an  in¬ 
trinsic  broadening  S/  due  to  screening  and  an  extrinsic 
alloy  broadening  S,,.  In  the  present  theory,  there  is  no 
need  to  suppose  that  the  exciton  suddenly  changes  its 
character  from  an  extended  effective-mass-like  state  to 
a  deep  localized  slate. 

We  need  to  calculate  Ej,  and  1=  S  /  +  asa  func¬ 
tion  of  alloy  concentration  v  The  calculations  are 
based  on  a  quantitative  coherent-potcntial-.ipprox- 
imaiion  (CPA)  hand  structure  Details  of  the  CP.A 
calculations  will  be  presented  elsewhere  Below,  we 
briefly  discuss  a  Green's-function  method  for  calculat¬ 
ing  Ej,  and  A  4 

The  one-particle  effective  Green's  function  in  CPA 
takes  the  form 

f;(E)  =  (E  /7  1(E)1  (II 

where  //  is  the  virtu. il-crysi.il  .ipproximaiion  Hamil¬ 


tonian  and  1(E)  is  the  self-energy.  The  site-diagonal 
Green’s  function  is  denoted  as 

FJE)^(,f>jG(E)\,i>a).  (2) 

where  ij)„  is  a  localized  orbital  of  specified  symmetry. 
Here  we  only  need  to  consider  a  =  s  for  Ai  symmetry. 
The  corresponding  function  in  pure  Si  is  denoted  as 
E®(£)  Following  the  theoretical  treatment  of  deep 
substitutional-impurity  levels,*  the  core-exciton  level 
for  pure  Si  is  determined  by 

E®(E)  =  (  T- (3) 

where  £,*'  is  the  site  potential  seen  by  an  5  electron  in 
bulk  silicon,  and  F  is  a  central-potential  parameter. 
For  a  chosen  value  of  Eq.  (3)  can  be  solved  for  £. 
and  vice  versa.  Then  £®“£°-£,  where  £®  is  the 
conduction-band  edge  in  pure  silicon.  Because  of  the 
uncertainties  in  the  value  of  experimental  £j®  and 
theoretical  F,  we  treat  E^  (or  F)  as  a  parameter.  The 
binding  energy  £»  m  a  Si;,Gei_;,  alloy  can  be  calculat¬ 
ed  by  solving 

£,(£)- 1  f  -  E,  1. (£))-'.  (4) 

where 

£,  =  x£/'+ (1 -x)E,"'.  (5) 

Then  £*  is  given  by 

E(,=  £r-E.  (6) 

The  calculated  values  of  the  conduction-band  edge 

and  the  exciton  level  measured  relative  to  the  top  of 
the  valence  band  are  plotted  in  Fig  I  The  band  gap 
increases  with  x  with  a  slope  discontinuity  at  =0.15 
The  dashed  lines  a.  h.  and  c  represent  exciton  levels 
oblaint  ’  with  E^®  =  0  1,  0  15,  and  0  .30  eV,  respective¬ 
ly.  The  binding  energy  E(,  is  also  an  increasing  func¬ 
tion  of  x,  with  a  slope  discontinuity  near  v  -0  15 
The  CPA  introduces  a  slight  bowing  m  E,  and  E^ 
Stnnati'®  has  c.ilculated  the  variation  of  A,  with  Ej, 
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FIG.  1.  Variation  of  the  band  gap  (solid  line)  and  the  Si 
2p  core-exciton  level  with  x  in  Si.Gei-,  alloys.  The  energy 
IS  measured  from  the  top  of  the  valence  band  The  dashed 
curves  a.  b.  and  c  repre.sent  exciton  levels  calculated  with 
T*  =0  1 .0.15,  and  0  3  eV,  respectively 


by  replacing  the  short-range  Coulomb  potential  with  a 
spherical  square  well  of  variable  depth  and  a  screened 
Coulomb  tail.  Strinati’s  results  can  be  used  to  estimate 
A,  corresponding  to  the  calculated  f*.  A/  decreases 
rapidly  with  £'*,  then  saturates  for  larger  £*. 

The  contribution  to  the  natural  linewidth  from  the 
alloy  broadening  is  calculated  by  a  consideration  of  the 
electron  pan  of  the  excilon  wave  function,  ij/.  The 
is  expanded  in  a  linear  combination  of  the  .spart  of  the 
conduction-band  wave  functions  <5 '. 

fl.S 

We  found  that  alloy  scattering  is  only  moderate  and  s 
scattering  is  dominant,  thus,  the  alloy  broadening 
A^l  £,')  IS  well  approximated  by 

A^(£)  =  .v(l  -  .v)fi,^Im£,(£).  (8) 

where  S,  is  the  difference  between  and 
Hence,  the  alloy-broadening  contribution  to  A  is  relat¬ 
ed  [0  the  alloy  broadening  of  the  band  states. 
AJk.E) 

A,  -  ~ 

^  k 

k  ft 

=  j’p,(£)A,  (£)(/£■ 

=  .x  (l  -  x  )f>;n  jt,;i  nj/:  (9) 

The  integral  in  Eq  (9)  is  ev.iluated  numericalK 


FICi  2  Variaiion  of  A  (solid  lines)  and  A/  (dashed  lines) 
with  X  for  three  values 


The  calculated  A.  which  is  the  sum  of  A^  and  A<.  is 
plotted  against  x  in  Fig.  2  for  three  values  of  E^.  In  all 
three  panels,  the  dashed  curve  represents  A,  and  the 
solid  line  represents  A.  It  is  seen  from  Fig.  1  that  the 
exciton  level  follows  the  A'  edge  of  the  conduction 
band.  Hence  the  binding  energy  relative  to  the 
conduction  band  edge,  remains  almost  constant  (for  a 
given  fj)  until  the  minimum  switches  from  the  \ 
edge  to  the  /.  edge  Because  of  the  change  in  the  slope 
of  Eg.  El,  decreases  rapidly  when  L  becomes  the 
minimum  Correspondingly,  A/  varies  slowly  until  the 
•V  to  L  crossover  and  then  increases  rapidly  This 
feature  is  clearly  seen  in  Fig  2 

For  Eft  =0  15,  the  A/  and  A.,  are  comparable  near 
v  =  0.50.  and  A/  dominates  for  all  small  \  and  large  v 
These  two  competing  mechanisms  give  a  relative 
minimum  near  .v==0.15,  a  broader  maximum  near 
A  =  0  50.  and  a  smaller  minimum  for  pure  silicon  As 
is  decreased,  the  relative  minimum  is  shifted  lo 
larger  x,  eg.  the  minimum  shifts  to  \  =  0  20  for 
/:/'=0  10  e\'  For  £(,‘’  =  0.15  eV.  the  positum  of  the 
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relative  minimum  is  in  agreement  with  the  experi¬ 
ment  (By  measuring  the  relatise  width  at  .v  =  0.15 
to  that  v  =  1.  one  can  make  a  better  estimate  of  F^'.) 
To  correlate  the  theory  with  experiment,  the  calculated 
1/A^  IS  compared  with  the  measured'^ 

in  Fig.  3,  where  is  the  edge  step  and 
( (i^ll  are  the  maximum  values  of  the  derivative 

of  absorption  spectra  with  respect  to  photon  energy 
Because  the  experimental  values  are  given  in  arbitrary 
units,  the  values  are  normalized  to  agree  at  .v  =  0.5 
The  observed  anomalous  behavior  near  .v  =  0. 15  and 
the  qualitative  x  dependence  in  that  region  is  clearly 
replicated  by  the  theory.  However,  the  calculation 
predicts  a  larger  maximum  at  a  =  1.  It  would  be  in¬ 
teresting  to  have  experiments  that  cover  the  entire 
range  of  x  to  further  test  thts  predictton 

for  larger  values  of  £j,  the  calculated  f(,  is  also 
large  and  hence  A/  decreases  slowly  with  x  Because 
the  broadening  is  determined  mainly  by  A.^,  the 
linewidth  is  expected  to  be  small  for  x  =0  and  x  1 
only,  this  occurs  for  £(['=  0  .K)  eV  For  negative  values 
of  /-.'jj*,  £(,  remains  negative  for  all  values  of  a  .Accord¬ 
ingly,  the  linewidth  is  broad  for  all  x,  and  there  would 
be  no  such  anomaly  as  in  Fig  3 
The  calculations  presented  m  this  Letter  are  slightly 
different  f  rom  alchemy  approximations  *  We  treat  the 
cer.'tral-cell  potential  Las  a  parameter  and  narrow  ns 
range  from  other  considerations  We  examine  values 
of  I  -  F,'"  of  -8  49.  -7  09.  and  -b  56  eV,  corre¬ 
sponding  to  Fl!  values  of  0  30.  0  15,  and  0,10  eV, 
respecticeK  If  the  strict  alchemy  approximation  were 
taken,  the  value  of  L-  F^'  would  be  F[-  -4  59 

e\'  in  the  tight-binding  approximation,  and  a  negative 


FK)  1  he  calculated  1  4'  values  (solid  line  I  compared 
with  ihc  cvperimenial  results  (marks)  from  Ref  17  The 
v.ilue  ol  calculated  4(Sii  is  0  17’’  eV  I  he  experiment  is 
lUirm.ili/'cd  In  (he  Iheorv  ,ii  i  -  0  s 


£®  ( - 0.10  eV)  would  be  obtained  When  long- 

range  interactions  are  included,  however,  the  above 
resonance  state  becomes  a  shallow  donor  level,  which 
IS  the  experimental  situation  for  a  P  impurity  in  Si 
Our  results  suggest  that  L  for  core  excitons  in  Si  is 
deeper  than  those  implied  by  alchemy  approximations 
However,  if  we  use  the  alchemy  approximation  as  a 
means  of  scaling,  the  value  of  L  for  Ge  3/r  core  exci¬ 
tons  should  be  deeper  than  that  for  Si  2/)  core  excitons 
Hence,  the  curve  corresponding  to  £/’  =  0.30  m  Fig  2 
should  be  a  reasonable  estimate  for  Ge  3pcore-exciton 
binding  energy  in  alloys.  Therefore,  we  do  not  expect 
to  see  an  anomalous  behavior  of  A  in  alloys  for  this 
case 

In  summary,  the  present  calculations  of  the  Si  2p 
core-exciton  binding  energy  and  linewidth  suggest  that 
the  exciton  level  is  about  0.15  ±0  05  eV  below  the 
conduction-band  edge  for  pure  Si  It  follows  the  X 
edge  for  a  >  0  15  in  the  Si;,Ge|_,  alloys,  and  £*  may 
eventually  reach  zero  in  the  dilute  limit  a  —  0  Our 
value  for  F^  represents  the  lower  end  of  the  previous 
measurements,'*''  but  is  in  good  agreement  with  a  re¬ 
cent  experimental'-*  value  of  0  120  ±0  03  eV  By  con¬ 
sidering  the  intrinsic  linewidth  and  the  alloy  broaden¬ 
ing.  we  can  explain  the  observed  relative  minimum  in 
the  linewidth  near  a  =0  15.  wuht'ut  requiring  a  sud¬ 
den  change  of  the  exciton  character.  On  the  basis  of 
this  calculation,  we  expect  the  corresponding  width  in 
pure  Si  to  be  even  smaller  than  that  observed  near 
A  =  0  15  We  further  argue  that  the  binding  energy  of 
Ge  3pcore  excitons  should  be  larger  than  that  of  Si  2p 
core  excitons  and  there  should  be  no  anomaly  in  the 
Ge  3p  linewidth  in  alloys 
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A  formula  for  alloy-scattering-limited  electron  mobility  in  semiconductors  is  obtained  for 
indirect  gap  systems  with  multiple  band  minima.  All  the  input  parameters  needed  are  defined 
explicitly.  The  drift  mobility  of  Si,  Ge,  , ,  which  has  a  dip  at  or  —0.1 3  and  a  broader  minimum  at 
X  —  0.5  is  calculated  by  adding  alloy  scattering  to  other  scattering  mechanisms  and  correlates  well 
with  the  measured  Hall  mobility. 


The  electron  and  hole  mobilities  in  semiconductors  are 
determined  by  the  band  structure  and  various  scattering 
mechanisms,  predominately  impurity  and  phonon  scatter¬ 
ing.  For  alloys,  the  mobility  is  also  affected  by  disorder  aris¬ 
ing  from  aperiodic  atomic  potentials  and  atomic  positions. 
Many  years  ago,  Nordheim'  and  Brooks*  obtained  an 
expression  for  alloy-scattering-limited  electron  mobilities  in 
metals  and  semiconductors,  respectively.  Brooks’  well- 
known  formula  reads 


42^ch*No _ 1 

3x(l  [AEffia'’ 


(1) 


where  jV„  is  the  number  of  atoms  per  unit  volume,  m*  is  a 
band-edge  effective  mass,  x  is  the  fractional  concentration  of 
one  of  the  species,  and  is  an  energy  parameter  character¬ 
izing  the  alloy  potential  fluctuations.  Although  this  formula 
has  been  widely  and,  to  some  extent,  successfully  used  for 
direct  gap  materials,  ’'’  the  identification  of  the  alloy  disor¬ 
der  parameter  AE  remains  uncertain.  Various  suggestions 
have  previously  been  made  for  AE.  e  g  ,  and  band-edge  dis¬ 
continuity"  or  band-gap  differences.’  Any  of  these  simple 
choices  is  bound  to  fail  when  one  applies  Eq.  (1)  to  more 
complicated  indirect  gap  systems  such  as  Si,  Ge,  ,  alloys, 
where  one  encounters  conduction-band  minima  transferring 
between  the  X  and  L  points  of  the  Brillouin  zone.  For  exam¬ 
ple,  if  is  taken  to  be  the  difference  in  corresponding  band 
edges,  then  one  finds  that  AE  —0.1  eV  for  the  X  iJ  )  valley 
and  1.2  eV  for  the  L  valley.  The  values  that  fit  the  experi¬ 
ment  are  about  half  this  value  for  L  and  -0.5  eV  for  T."  The 
purpose  of  ihis  letter  is  to  resolve  the  identity  of  for  indi¬ 
rect  gap  materials. 

Moreover,  there  is  a  problem  with  the  rrt*  that  enters 
Eq.  ( 1 1  For  direct  gap  alloys,  the  band-edge  etfective  mass  ai 
r  naturally  enters  Eq.  ill.  For  the  indirect  gap  alloys,  the 
effective  mass  is  anisotropic  and  hence  an  appropriate  mass 
must  be  chosen.  Previous  authors"  ’  have  chosen  m*  to  be 
the  effective  conductivity  mass  m*.  We  shall  show  that  dif¬ 
ferent  masses  enter  for  different  cases. 

The  first  unambiguous  assignment  for  J  fin  a  direct  gap 
alloy  was  given  by  Hass  et  a!.'  1  o  estimate  the  limiting  elec¬ 
tron  mobility  in  Hg,  Cd,  ,  Te  based  on  a  tight-binding  iTBi 
band  description,  they  defined  AE  to  be  fAE..  where /,  is 
the  s  fraction  in  the  density  of  stales  and  Jf,  is  the  difference 
between  the  s  atomic  term  values  of  the  Hg  and  Cd  atoms.  Bv 
extending  this  approach  to  alloys  with  indirect  gaps  and 


multiple  bands,  we  show  that  all  the  uncertainties  identified 
above  are  resolved.  Our  generalized  Brooks'  formula  will 
then  be  applied  to  Si.  Ge,  _ .  systems  to  explain  their  ob¬ 
served  mobility.’’* 

Because  Brooks'  formula  has  never  been  derived  expli¬ 
citly  in  the  literature,  we  rederive  it  first  and  then  generalize 
it.  Consider  the  case  of  a  single  band  with  an  isotropic  effec¬ 
tive  mass.  The  dc  electronic  conductivity  based  on  the  linear 
response  theory’  is  given  by 

where  the  energy-dependent  a\()  in  the  weak  alloy  scattering 
limit  is 


<7(f)  =  (eV3)i’V)Z)(f)Tte). 


(3) 


D  (f)  is  the  density  of  states  (DOS)  per  unit  volume  for  both 
spins,  so  D  (f )  =  with  —  N^/l  being  the  number 

of  unit  cells  per  volume  (for  the  diamond  structure,  half  the 
number  of  atoms  Nc,  per  unit  volume)  and^)f  I  being  the  DOS 
per  unit  cell  per  spin.  The  mean  square  velocity  i''(f)  for 
carriers  with  energy  e  is  given  by 


r  pie) 

The  scattering  lifetime  for  carriers  with  energy  f,  riel,  is  re¬ 
lated  to  the  alloy  broadening  A  If)  by  rtf  I  =  ft/lA  if  i.  where 
the  energy  J  (f  |  is  the  imaginary  part  of  the  self-energy  in  the 


averaeed  allov  Green’s  function.  For  w  eak  scattenne  J 


J  If)  =  17x11  —  xl  \AE  r’^lfi. 


I5l 


where  in  a  tight  binding  iTBi  description  AE  is  the  difference 
in  the  term  values  of  the  constituents.  Then  the  mobility  is 
p,  =  a/ne  with  the  electron  density  given  by 

n  =  2  J/ifg7|f  )(ff.  lb) 

For  a  nondegenerate  semiconductor. /if  i  is  ihe  Boltzmann 
distribution  and/lf  I ->  e  '  Furthermore,  for  a  para¬ 

bolic  bandfiki  =  fi'k  V2m*,/)|f|  =  i2m*i'  / '  ’/4r"'^-',then 
all  the  above  equations  can  be  combined  to  arrive  at  Eq  1 1 1 

For  a  real  semiconductor  alloy  in  a  TB  description,  the 
alloy  scattering  can  be  characterized  by  two  parameters 
and  AE^..  the  differences  in  s  and p  term  values  between  two 
subslitutional  atoms  Then  an  efTective  broadening  is  given 
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1  ABLE  I  Calculation  paramciers. 


d  (f)  =  (J,^,  +  (7) 

where/),,/?,  are  partial  density  of  states  (PDOS)  and  J,  and 
J,  are  similar  to  Eq.  |^),  with/?  replaced  by  p,  and p,.  re¬ 
spectively.  For  Hg,  ,Cd,Te,  the  s  disorder  is  predomin- 
ent’ *"  and  one  can  neglect  4£,.  Defining  p,  = /,p  (and 
p,  =/fP],  one  arrives  at 

id=:)Tjc(l  -x){f,^E,fp. 

Thus  as  was  pointed  out  by  Hass  et  al.'f.dE,  plays  the  role 
of  AE  in  this  special  case  where  AE^  can  be  neglected. 

For  an  alloy  with  a  single  indirect  gap  minimum,  one 
has  to  consider  both  s  and  p  contributions  to  the  alloy  broad¬ 
ening  and  the  masses  that  enter  p  and  o'.  Again,  Eqs.  (2Hfi| 
can  be  combined  to  yield 


_ _ _ 

[3x(l  -x]n,*m:(mr]''\kT]''^\J/lAEl+/lAEl)]' 

(8) 

where  m*  and  m*  are  respectively  the  longitudinal  and  the 
transverse  mass  at  the  band  edge,  and  A^,.  is  the  number  of 
equivalent  minima,  e.g.,  6  for  Si.  The  conductivity  mass  m* 
comes  from  averaging  y*  in  Eq.  (2)  and  is  given  by  3(2/»i* 
+  l/m^l  '.  Equation  (8)  clearly  identifies  the  masses  and 
the  energy  parameter  that  enter  Brooks'  formula. 

Next  we  consider  a  still  more  complicated  case  where 
the  contribution  to  the  mobility  comes  from  more  than  one 
band.  For  example,  in  Si,  Ge,  _ ,  the  -IT  and  L  minima  cross 
near  jr  =  0. 1 5. "  There  are  now  two  contributions  to  the  net 
conductivity,  so  =  2<t.,  where  /  is  X  or  L.  The  quantities 
w,V).  D,  (€),  and  N,  (f )  now  take  different  values  for  different 
bands.  The  structure  of  r,  (e)  requires  more  careful  consider¬ 
ation.  The  complication  comes  from  the  fact  that  the  effec¬ 
tive  broadening  A  is  still  given  by  Eq.  (7),  but  p.,p,.  andp 
contain  contributions  from  both  the  bands.  The  proper  ex¬ 
pressions  arep  =  2,p,A'|.  andp„  =  where  /  =  X 

orL.o  =  rorp.  and.V,*^  =  6,  A'f-  =  4.  The  equation  for  J  is 

A  (f|  =  7rx{]  -  x}^(^^/,„y:p,{f)AE.,  jj^2^p,(e)A’:J.|9| 

The  mobility  associated  with  the  /th  band  is  defined  as 
p,  =  a,/{n,e).  then 

p,  = - - /, .  ilOl 

3x(  1  -  .rl  -  ], 


Parameter 


fp,if)(y.v:p^(fi\'  •  *' 

- 


”  df.  iiii 


Thus,  the  generalized  formula  no  longer  has  the  explicit 
X  and  r  dependences  of  the  original  Brixtks'  form.  Howe'er, 
all  the  quantities  needed — the  masses,  the  scattering  param¬ 
eters  J  A',,  .  the  band  gaps,  and  the  fractions —can  be  e'al- 
uated  iheoretically  without  resorting  to  experimentally  fit¬ 
ted  parameters.  To  demonstrate,  we  shall  apply  Eq.  (lOl  to 
Si,  Ge,  . .  The  band  quantities  are  obtained  from  our  CPA 
calculation.  “We  found  that  the  effective  masses  varv  weak- 


Si.Ge,  ,  sysicms 
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ly  with  the  concentration,  so  and  m*  are  assumed  to  be 
constant  and  assigned  the  values  0.97  and  0.19  for  the  X 
minima  and  1.64  and  0.082  for  the  L  minima,  respectively. 
The  calculated  energy  gaps  for  the  X  [A  |  follows  the  func¬ 
tional  form  E'*'  =  a  +  bx  +  ex'  and  for  L  is  given  by 
=  A  +  Bx  +  Cx^.  All  the  parameters  of  our  calculations  are 
listed  in  Table  I. 

To  correlate  the  calculation  with  the  measured  mobili¬ 
ties,  we  need  to  have  an  estimate  of  scattering  rates  l/r„  due 
to  impurities  and  phonons.  A  crude  approximation  is  to  as¬ 
sume  l/r,,  for  a  given  valley  to  be  the  same  as  the  appropriate 
constituent’s  values  and  add  to  it  the  alloy  scattering  rate  1/ 
.  Then  the  average  mobility  and  the  mobility  from  the  ith 
minimum  in  the  alloy  are 

=  +  (12) 
p*  is  given  by  Eq.  |10)  and  p°  are  the  measured  drift  mobili¬ 
ties  for  Si  or  Ge.'^  The  drift  mobility,  calculated  from  Eq. 
( 1 2),  is  plotted  as  a  function  of  alloy  concentration  x  in  Fig.  1 . 

For  x<0.05  and  .x>0.20,  the  energy  difference  between 
the  X  and  L  edges  is  large  enough  so  there  is  a  negligible 
contribution  to  the  mobility  from  the  higher  minima.  In  the 
Si,  Ge,  _ ,  system,  the  5  scattering  is  predominent.  Because 
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the  L  edges  have  a  larger  s  content,  alloy  scattering  reduces 
the  average  mobility  substantially  for  small  x.  Even  though 
the  i  content  is  almost  the  same  for  all  x  >0.20  at  the  Jlf  edge, 
the  mobility  still  decreases  to  0.5  as  shown  in  Fig.  1  be¬ 
cause  of  the  jt(  1  ~  x)  term  in  Eq.  ( 10). 

An  interesting  feature  is  obtained  for  the  compositions 
O.Bcjc  <0.18.  The  average  mobility  attains  a  local  mini¬ 
mum  near  x  =  0.14  and  a  smaller  maximum  at  x  =  0.17. 
This  feature  occurs  because  of  the  X  to  L  crossover."  For 
x<0. 14,  the  major  contribution  to  fx  comes  from  L  minima. 
Near  x  =  0. 14,  the  density  of  states  increases  because  the  X 
and  L  minima  merge.  So  the  alloy  scattering  increases  there 
and  the  average  mobility  decreases.  Forx>0. 14,  the^bands 
have  the  lower  minima.  As  the  s  content  is  small  at  the  X 
minima,  the  reduced  alloy  scattering  increases  the  average 
mobility.  For  larger  values  of  x,  the  jr(  1  —  x)  term  takes  over 
and  the  mobility  varies  as  shown.  The  values  of  measured 
Hall  mobility  in  Si,  Ge,  ,  systems  are  also  plotted  in  Fig.  1 . 
The  interesting  feature  near  x  =  0. 14  is  clearly  seen.  Since 
the  experimental  drift  mobility  //„  for  Si,Ge,  _  ,  is  not 
available  and  the  generalization  of  Eq.  ( 10)  to  Hall  mobility 
fx„  is  less  clear,  we  present  the  calculated and  experimen¬ 
tal (Ref.  7,8)  here.  While  we  do  not  expect  quantitative 
agreement,  because can  range  from  1  to!,""  we  do 
expect  them  to  display  the  same  qualitative  x  dependence.  It 
is  rewarding  to  note  the  similarity  in  the  trend  in  Fig.  I. 
Previous  authors  explained  the  dip  in  the  mobility  curve  by 
including  intervalley  scattering  with  an  arbitrary  adjustable 
coupling  constant.'’  Our  calculations  automatically  include 
that  portion  of  intervalley  scattering  that  results  from  alloy 
disorder  with  a  coupling  constant  set  by  the  atomic  proper¬ 


ties  of  the  constituents  However,  the  additon  of  intervalley 
scattering  mediated  by  phonons  and  impurities  is  expected 
to  increase  the  dip  near  crossover. 

In  summary,  an  expression  for  alloy-scattering-limited 
charge  carrier  mobilities  is  derived  for  indirect  gap  alloys 
with  multiple  bands.  This  expression  reduces  to  Brooks'  for¬ 
mula  for  direct  gap  alloys.  The  quantities  m*  and  AE  can  be 
calculated  exactly.  Alloy  scattering  accounts  for  the  ob¬ 
served  mobility  features  in  the  Si.Ge,  alloy,  including 
the  anomaly  near  the  L  to  AT  (d  )  crossover. 
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ABSTRACT 

Several  recent  theoretical  studies  of  the  local  structure  of  semiconductor  alloys  are 
summarized.  First,  dilute  limit  calculations  of  local  bond  lengths  and  mixing  enthalpies 
are  discussed.  These  calculations  include  effects  due  to  both  bond  length  and  bond- 
angle  distortions,  as  well  as  local  chemical  rearrangements.  Then,  a  new  statistical 
theory  of  concentrated  alloys  is  described.  Deviations  from  random  alloy  distributions 
(microcluslers)  are  predicted. 


INTRODUCTION 

This  paper  summarizes  our  recent  theoretical  studies  directed  toward  understand¬ 
ing  the  microscopic  structures  of  pseudo-binary  semiconductor  alloys  A,B,.,C.  We  first 
present  a  detailed  calculation  of  the  local  bond  length  relaxation  in  the  dilute  limit 
X  -•  0,  i.e.  the  case  where  an  .A  atom  is  substituted  for  a  B  atom  in  a  BC  compound. 
The  mixing  enthalpy  parameter  fl  is  found  to  be  related  to  small  excess  substitution 
energies.  These  excess  energies  are  calculated  directly  through  a  minimization  pro¬ 
cedure.  Thus,  the  accuracy  of  the  predicted  0  is  not  limited  by  trying  to  find  small 
differences  between  large  numbers.  The  theory  is  then  generalized  to  concentrated 
alloys  using  statistics  based  on  combinations  of  tetrahedral  clusters  of  five  atoms.  Our 
results  predict  that  microclustering  occurs  in  a  majority  of  alloys.  We  conclude  by  iden¬ 
tifying  systematic  correlations  between  the  theory  and  several  experiments. 

Before  discussing  the  calculation,  it  is  useful  to  provide  some  background  about 
the  structure  of  these  alloys.  It  was  customary  to  assume  that  these  alloys  have  two 
sublatlices  in  which  the  C  atoms  occupy  one  sublattice,  and  and  B  atoms  are  ran¬ 
domly  distributed  on  the  other.  This  picture,  referred  to  as  the  virtual  crystal  apprc'Xi- 
malion  IN'C.A),  implies  that  the  nearest-neighbor  (nn)  bond  lengths  in  the  alloy  are  the 
concentration  weighed  average  values,  i.e.  d,(-  =  dgo  =  d  =  x  dj^^  -t-  (1-xl  d^o  where 
the  values  with  a  superscript  (01  denote  the  pure-crystal  values.  On  the  other  hand, 
according  to  Pauling's  covalent  r.adii  approximation  (CR.\),  the  local  bond  lengths 
retain  their  respective  pure-crystal  values,  i.e.  d^c  =  d|^  and  dg,-  =  d^®' 

If  we  define  f'g  =  (d  -  d^^')  /  d  and  b  —  (<i  -  d^c)  /  f*-  'ben  the  ratio  bjb^  in  \'C \ 
is  rero,  hilt  in  ('R,\  it  is  I.  However.  Mikkelsen  and  Boyce'"  found  from  their 
evpcrini'  nt  ivii  <  .a.,lni..^  Vs  that  the  nn  bond  lengths  do  not  fit  either  \  (  '.-\  or  t  'l'  A. 
Instead,  they  found  the  value  of  blt'^  to  he  close  to  3/4.  Since  then,  similar  experime-.i-- 
have  been  done  for  a  number  of  rinc-Men<ie  p.seudc>-bin.ar.v  alloys.  ’•  an'l  the  ,3,  t  -uh 
a(3*ears  to  t>e  ‘(uite  general. 
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DILUTE  LIMIT 

The  dilute  limit  is  the  easiest  case  but  is  still  not  trivial.  Its  solution  provides 
both  end-point  results  (x  =  0  and  I),  as  well  as  insight  into  the  exteasion  to  the  concen¬ 
trated  alloy  C3M\  A  complete  description  of  this  case  is  being  reported  elsewhere;  here 
we  summarize  the  essential  results.  The  substitution  energy  A,  for  an  A  atom  replac  ng 
a  B  atom  in  a  BC  compound  is  calculated  and  minimized  to  find  the  relaxed 
configuration.  A,  can  be  written  as  A,  =  4  (AE^  +  AE,  -t-  AE^^),  where  AE^  is  the 
binding  energy  difference  between  the  AC  and  BC  compounds.  AE,  is  the  strain  energy, 
and  AE,|,  is  a  chemical  energy  shift.  All  AE'  s  are  energies  per  bond.  Then. 
AE  =  AE,  +  AEch  is  the  excess  energy  per  bond  for  the  impurity  substitution.  AK,  is 
cale\ilated  by  dividing  the  crystal  into  two  regions.  Outside  R  (which  is  the  distance  of 
the  second-shell  atoms  to  the  impurity),  the  distorted  crystal  is  treated  as  an  elastic  con 
tinuum  with  a  radial  displacement  field  which  is  inversely  proportional  to  the  square  of 
the  radial  distance,  so  AE,*®"''*  =  l/4RCu’,  where  C  is  an  effective  shear  coefficient. 

C  =  «  (1.6  (C„  -  C,j)  -H.S  C,,),  (1) 

and  u  is  the  magnitude  of  the  displacement  at  R.  Inside  R,  the  strain  energy  A,*'"’  is 
treated  with  a  valence  force  field  (\TF).***  Finally,  the  chemical  energy  shift  AF.js  is  cal¬ 
culated  from  Harrison's  model  and  arises  from  changes  in  the  metalization  energies 

caused  by  different  bond  lengths  Ad  =  dgc  -  d/,c  covalent  energies 

AVj  =  Vj(AC)  -  V3(DC).  Note  that  =  (d^  -  d|®ci)  /  d^'  and 

^  •  dj^(^)  /  d^°'  in  this  dilute  limit,  so  the  excess  energy  AE  can  be  expande.l  up 

lo  M’cond  order  in  6.  u,  and  AV3  For  a  given  pair  A  and  B,  is  an  explicit  furiciit>ii  f'f 

6  and  u.  Minimization  of  AE  with  respect  to  6  and  U  leads  lo  the  equilihriujn  lt>ral 
bond  length  d^c  energy  AE.  Then.  AE  is  used  to  estimate  the  mixing  enthalpy 
parameter  0  in  the  mixing  enthalpy  =*  xfl-xjO  b\ 

n  =  2  (AE  (A  in  BC)  +  AE  (B  in  .XC)).  (21 


A  syMematic  comparison  with  other  models  based  on  strain  ♦•nergy  alone  shov^v 
that  an  incr»*a.se  of  the  range  of  the  fixed  boundary  R  increases  ifif  relaxation  of  d^(  .  i.e. 
It  causes  6/6q  to  incrpa.se.  The  iiichisjon  of  the  bond  angle  restoring  force,  on  the  ot  her 
han<l.  reduces  the  relaxation.  It  turns  otjt  that  a  delicate  canceilalion  of  these  tv^f^ 
etTerl.s  causes  3  simple  spring  model  j»oinled  out  by  Shih  el  ai.  (SsHS)'''  l«'  vieM  accu. 
rate  results.  In  this  motlej  *1  'xhere  .)  anti  (*|  are  ihe  l.'.tui- 

stretching  force  constants  for  the  h<xst  (BC)  and  the  imf'uriiy  lAC^  crystals, 
ci  ^  a|.  this  model  predicts  (>  --  for  a  zinc-blende  alloy.  .XJihough  our  full  pertur- 

hat  ion  l  hoory  '  Ff'T)  and  the  \'F'I-  nioiiej  (,f  Nh»riin  and  hunger  i  M/  predin  t  with 

an  averaee  abst-iule  deviation  comparable  t<>  the  experinienta!  uricertainty  ryf  itii 

'■iinpie  sprim:  m  del  e\#*n  better. 

.’h  it  V*  fiih'  >h  'Mir  ih'-’  r>  .  N?/  and  ''SJIS.  th'-  va.:ic'  rire  directly  -'  if, 
ialcfj  tji  :i;;y  a-j]us»:ibl«  jaraiin’ers.  c*ur  theory  and  -iCf*''  \Mlh  the  r 

niei’t  .L‘-  w*!l  •  -■  e\ri;  ilightly  ’h.an  the  one-parameter  (he^  nes. .Mth'-ugh  ")'■ 

theory  pre-lie'*-  a  i.egativi  U  vuiue  f'‘r  ail  three  ;(ia.  W)  ali'Vs.  the  macnilui*- 
fi.17  kc.'ij/im  it  1  is  t*).^  small  m  at  c-.ium  for  the  .-rdering  v-f  grown  at 

fi(X)  to  7(X)’('  found  recently.  The  caJciiiated  U  values  a!s<y  provide  guidance  in 
separating  miscii«le  fn-m  iinml-cibU-  ailoys.  In  a  random  alloy,  ihe  criterion  for  alloy 


i 
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mixing  for  all  x  is  T  >  T„  where  Ibe  critical  temperature  T  is  given  hy  T,  =  n/21<, 
with  R,  being  the  universal  gas  constant.  Figure  1  is  a  plot  of  TJTj  against 
1  *0  I  /  I  *«i  I  '  Tj  is  the  lower  of  the  two  constituent's  melting  temperatures,  and 
4^  =  1.63Xib  '•'ifh  Xo>  being  the  ratio  of  rms  bond  length  amplitude  fluctuation  to  the 
bond  length  at  Tj.  The  simple  spring  model  gives  T^/Tj  =  indicated  by  the 

solid  curves. 

There  is  an  empirical  rule  stating  that  a  miscibility  gap  will  occur  if  |  Ag  I 
between  two  alloy  components  exceeds  7.5^.  However,  if  T^/T^  is  plotted  against 
1  A)  I  .  the  simple  spring  model  would  not  exhibit  a  smooth  simple  quadratic  curve,  and 
our  theoretical  points  would  be  much  more  scattered.  This  suggests  that 
1  ^0  I  /  I  1  >  1  is  a  better  criterion  than  |  4  |  >  0.075.  Figure  2  also  clearly  shows 
the  chemical  effects,  namely  negative  and  positive  chemical  energies  AEj,,  for  cation  and 
anion  substitutions  respectively.  The  full  theory  and  the  experiments  correlate  within 
the  experimental  uncertainties.  The  simple  SSHS  model  clearly  is  an  excellent  universal 
representation.  However,  T^/T,  varies  faster  than  quadratically  for  larger  |5o/^ml 
values,  as  born  out  from  both  the  experimental  data  and  the  full  theory. 


CO.NCENTRATED  ALLOYS 

Turn  now  to  the  concentrated  alloy  case.  First,  an  improved  statistical  model  is 
required.  We  have  extended  regular  solution  theory  based  on  pair  energies  to  one  for 
five-aiom  clu.sters.  For  an  alloy,  the  building  blocks  are  clusters  of 

A(m)n(-l-m)C,  where  m  ranges  from  0  to  •!.  For  a  given  alloy  concentration  x  and  for  a 
given  set  of  energies  <„  atwociatod  with  these  clusters,  we  have  derived  expressions  for 
the  cluster  population  distribution  x„  =  /  N.  where  is  the  total  number  of  unit 

cells  and  n,„  is  the  averaged  number  of  cells  with  A(m)B(‘l-m)C  clusters.  The  partition 
function  Z  is  obtained  using  a  steepest  descents  argument  which  then  yields  the  mixing 
llelmhott  free  energy  AF.  The  result  reduces  to  Guggenheim's  tetrahedron  case*''  if 
pair  potentials  (for  the  second-neighbors)  are  assumed.  .Another  major  difference  is  that 
w"  only  need  to  solve  a  single  quartic  equation,  while  Guggenheim  needed  to  solve  four 
>imultaneous  quartic  equations. 

The  key  to  the  problem,  however,  lies  in  the  calculation  of  the  energies  If  one 
assumes  that  the  site  of  the  leirahedra  for  all  m-clusters  at  a  given  alloy  concentration 
takes  on  the  corrc.sponding  \  C.A  values  but  allows  tbe  central  C  atom  to  relax,  then  the 
.'in  riles  .V.  functions  of  x  behave  like  those  shown  in  Figure  'i(a).  There  are  at  lea.st  two 
:i,,ijor  in  ibis  resuli.  First,  the  energies  ar-  t.vs  large  and  would  correspond  In  f! 

\aiues  mariv  limes  the  experimental  values.  .Second,  at  x  =  0.7.5,  0.5,  and  0.2.  tbe.se 
(‘nergies  iriiplv  conifs/und  formation  for  .A^BiC^.  .A-iB^C  4.  and  .^58304,  respectively, 
which  !«  opposite  to  the  known  tendency  for  splnodal  decomposition  of  Ga,Ini.,.As  at 
'.IW  I.  However,  if  the  local  cell  volume  c>f  each  cluster  is  allowed  to  be  in  mechanical 
•'quili!  ri'im  with  a  continuous  medium  with  an  effective  shear  coefficient 
r  -  :l  vM'Mi.  “hiTc  the  (.'  value  for  the  pure  material  is  given  hy  F.q,  III. 

then,  -he  corre--pondinR  energies  a.s  3  function  'if  x  are  given  in  Figure  -.'(bl,  whicli 
;i.  w  vields  a  re'v-onable  value  of  niixing  entbalpo  and  correctly  predicts  the  tendency 
■  ward  spinodal  decomposition  at  low  temperature.  With  this  set  of  energies,  one  can 
t.di  aiculale  the  c|,jst,.r  distribution  xm-  and  compare  them  with  the  corresponding 
values  for  a  random  alloy,  i.e,  x'"'  =  (d^tx'^d  x)*^"'.  Figure  li  shows  the  deviation 
fr-m  randomness  Ax„  =  x,„  -  x^®'  as  a  function  of  x  for  four  arbitrarily  chosen  growth 
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FIGURE  1  PLOT  OF  T^y  Tj  AS  A  FUNCTION  -  hf,  FOR  T^  OBTAINED  FROM  (>)  PRESENT  THEORY 
AND  (bl  THE  experimental  Q  IN  Ul.  THE  DOTS  ARE  FOR  ANION  ALLOYS  AND  CROSSES 
FOR  CATION  SUBSTITUTION  THE  SOLID  LINES  IN  BOTH  PARTS  CORRESPOND  TO  THE  11 
BASED  ON  The  SIMPLE  THEORY  OF  REF  3  THE  DASHED  LINES  AT  T^/ T2  •  OSEPARATE 
THE  MISCIBLE  FROM  IMMISCIBLE  GROUPS 


Figure  2  the  cluster  energies  as  a  function  of.  calculated 
IN  TWO  DIFFERENT  WAYS  AS  0»SCUSSED  IN  TmE  TEXT 
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ABSTRACT 


Starting  from  realistic  band  structures  of  the  constituent  materials,  the  electronic 
structure  of  Si^Gei.x  alloys  are  obtained  in  the  coherent  potential  approximation  (CPA) 
Various  quantities,  including  the  bowing  parameter  of  the  fundamental  gap  and  the 
energies  of  several  optical  gaps,  the  masses,  and  the  linewidths  of  the  Eg  and  E]  transi¬ 
tions,  are  calculated  on  the  basis  of  both  diagonal  and  oif-diagonal  CPA.  All  of  the 
band-energy  and  line-width  predictions  are  in  good  agreement  with  experiments.  Fur¬ 
thermore,  the  theory  yields  alloy-scattering-limited  electron-drift  mobility  in  qualitative 
agreement  with  experimental  results. 


I.  INTRODUCTION 


Semiconductor  alloys  offer  the  freedom  to  design  material  properties  by  choosing 
appropriate  alloy  constituents.  In  some  cases,  the  physical  properties  of  the  alloys  can 
be  quite  different  from  those  of  the  constituents.(^'^)  In  recent  years,  there  has  been  a 
renewed  interest  in  SixGei.j,  alloys(^)  and  superlattices.(^^®)  Because  silicon  is  the  most 
technologically  advanced  semiconductor,  the  results  on  SiOe  systems  have  many  poten¬ 
tial  applications. 

The  lattice  constants  of  silicon  and  germanium  differ  by  ~4%.  Hence,  the  strain 
introduced  in  the  formation  of  Si^Gej^^  alloys  can  affect  the  band  structurei^®)  and  the 
transport  properties. (^)  Prior  authors  used  virtual  crystal  approximation  (VCA)(*U12) 
and  coherent  potential  approximation  (CPA)(^^)  to  study  the  band  structure  and  related 
properties.  Either  because  of  less  accurate  band  structures  of  the  constituent  materials, 
or  because  of  the  approximations  involved  in  the  alloy  formalism,  these  calculations 
predicted  only  trends  of  specific  quantities,  not  quantitatively  accurate  results.  Because 
the  s-state  site  potentiab  (c,)  for  silicon  and  germanium  differ  by  approximately  1.5  eV, 
VGA  cannot  accurately  describe  effective  masses  and  other  finer  details  of  the  band 
structure.  Because  of  the  use  of  poor  basis  functions,  earlier  CPA  work(^^)  predicted 
alloy  broadening  of  conduction  band  states  substantially  differing  from  experiment.  The 
purpose  of  this  paper  is  to  correct  these  flaws  and  treat  transport  phenomena. 

Because  of  a  substantial  difference  between  the  site  potentiab  and  lattice  constants 
of  silicon  and  germanium,  we  incorporated  both  chemical  and  structural  disorder  in  the 
calculation  of  the  electronic  structure  of  Si^Gcj.,  alloys.  Thus,  both  diagonal  and  off- 
diagonal  CPA  are  included  in  the  predicted  band  structure  and  related  quantities.  Parts 
of  the  band  structure  have  been  used  to  study  the  Si-2p  core-excitonl^"^)  and  the  alloy 
mobilities. A  comprehensive  report  of  the  calculations  and  results  is  presented  here 
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The  rest  of  the  paper  is  arranged  as  follows.  The  detailed  procedure  of  htting  sili¬ 


con  and  germanium  band  structures  is  given  in  Section  II.  The  VCA,  CPA,  and  off- 
diagonal  CPA  calculations  are  described  in  Section  III.  The  results  and  interpretation  of 
the  alloy  band  structures  and  mobility  are  given  in  Section  IV. 


II.  BAND  STRUCTURE  BASIS 

In  order  to  derive  an  accurate  alloy  band  structure,  one  must  start  from  a  realistic 
band  structure  of  the  constituent  materials.  Chen  and  Sher  have  developed  a 
method(^^)  following  a  prescription  of  Kane(^^)  and  Chadi(^^)  that  includes  all  long- 
range  interactions,  and  then  fine  tuned  the  band  structure  with  an  adjustable  local  Ham¬ 
iltonian.  Because  the  details  have  already  been  published, the  underlying  method 
will  be  presented  here  in  brief. 

Gaussian  orbitals  of  the  type  a  (  a  can  be  s,  p^,  p^,  or  p,)  for  each  sublattice  in  a 
cell  are  used  to  construct  the  corresponding  Bloch  basis.  In  this  basis  set,  the  overlap 
matrix  and  the  Hamiltonian  derived  from  empirical  pseudopotentiak  can  be  calcu¬ 
lated. It  is  possible  to  cast  the  problem  in  a  basis  set  of  Gaussian  orbitals  in 
which,  in  crystal  units  (cu),  the  same  exponential  factors  apply  for  all  Ill-V  com¬ 
pounds. In  this  universal  basis,  the  overlap  matrix  and  the  kinetic  energy  matrix  are 
same  for  all  III-V  compounds.  Then,  by  a  unitary  transformation,  the  basis  set  is  ortho- 
normalized. (“*^)  The  Hamiltonian  in  this  new  basis  set  is  denoted  Ho(k).  The  band 
structure  resulting  from  this  method  reproduces  the  results  of  elaborate  band  structure 
calculations  within  a  few  percent  throughout  the  Brillouin  zone  (BZ),  To  establish  accu¬ 
rately  certain  important  band  structure  features  adjacent  to  the  gap,  an  extra  small  8x8 


Hamiltonian  matrix  Hi(lc)  is  added  to  Ho(k).  This  Hi(k)  has  the  form  of  a  tight-binding 
(TB)  Hamiltonian,  in  which  only  the  nearest  neighbor  interactions  are  included,  and 
stimulates  the  effect  of  nonlocal  pseudopotentials  and  an  expanded  orbital  set.  The  total 
Hamiltonian,  H(ir),  in  this  orthonormalized  basis  set  is  diagonalized  to  obtain  the  band 
energies  and  the  corresponding  wave  functions. 

Following  this  procedure  with  the  same  exponential  factor  0  —  0.26  in  the  Gaus¬ 
sian  orbitals  for  both  silicon  and  germanium,  the  matrix  Ho(k)  is  obtained.  For  silicon 
and  germanium,  H^  contains  6  adjustable  parameters:  namely,  the  corrections  to  the 
term  values  A,  and  Ap  and  to  the  nearest-neighbor  interactions  V„,  V.p,  v„,  and  Vj,y. 
The  values  of  A,  Ap^  V„  and  are  determined  from  fitting  the  three  experimental 
energy  gaps(21-28)  at  r(k  =  0):  Fj*  -  Fi,  Fjs  -  F2S'  ,  and  Fj/  -  Fjs*  ,  and  the  photo¬ 
electric  threshold  (PT)  values  -5.07  and  -4.80  eV  for  silicon  and  germanium  respec¬ 
tively  .(2^)  The  remaining  parameters  V,p  and  arc  obtained  from  the  experimental 
values(2 1*28,30)  of  gaps  X,,  -  X,y  and  Ljj  -  Ly  Some  adjustments  in  these  input 
quantities  are  made  to  obtain  an  overall  good  band  structure  with  more  accurate 
effective  masses.  Table  I  lists  the  empirical  pseudopotential  form  factors  and  the  param¬ 
eters  used  to  obtain  the  band  structure.  The  calculated  band  structures  and  experimen¬ 
tal  values  are  given  in  Table  II.  From  Table  II,  one  can  see  that  an  excellent  fit  to  the 
silicon  and  germanium  band  structure  is  obtained;  All  the  calculated  values  lie  within 
the  experimental  uncertainties.  The  optical  difference  between  Lj^  and  Ls,,  Fjj'  and  F^ 
are  in  excellent  agreement  with  the  known  optical  transition  values. 

Although  the  calculated  effective  transverse  masses  agree  very  well  with  experi¬ 
ment,  the  effective  longitudinal  mass  for  germanium  is  less  than  the  experimental  value. 
This  is  due  mainly  to  our  attempt  to  have  a  common  0  and  the  choice  of  local  pseudopo¬ 
tentials,  causing  HqO')  same  in  crystal  units  (cu)  for  both  germanium  and 


silicon.  Because  of  the  common  Ho(k),  the  alloy  disorder  is  contained  in  these  adjusted 
parameters.  This  Hq  would  also  be  useful  for  the  interface  and  superlatticel^O)  prob¬ 
lems.  If  we  grant  ourselves  the  freedom  to  adjust  longitudinal  effective  mass  in  ger¬ 
manium  can  be  fitted  to  the  experimental  value.  When  is  changed,  the  Lie,  W  ,  will 
also  change.  We  have  chosen  not  to  do  this  because  little  is  gained  for  the  extra  com¬ 
plexity.  For  an  indirect  gap  semiconductor,  the  important  effective  mass  used  in  tran¬ 
sport  studies  is  the  conductivity  mass,  3  (l/m*,|  +  2/m*ei)  ■  Because  m*,)  »  m’^t  in 

germanium,  m*e  will  not  be  much  different  if  a  less  accurate  value  of  m*,!  is  used.  More¬ 
over,  the  Si^Gei_„  alloys  which  have  potential  device  applications  are  in  the  silicon-rich 
region,  where  the  effective  mass  at  the  L  edge  is  not  expected  to  affect  the  further  stu¬ 
dies. 

It  is  important  to  note  that  an  excellent  fit  to  the  experimental  values  can  be 
obtained  with  only  seven  adjustable  parameters  {0,  A„  Ap,  V„,  V,p,  V„,  V^).  with  0 
being  universal  in  cu.  The  calculated  band  structure  of  silicon  and  germanium  are 
shown  in  Fig.  1(a)  and  1(b)  respectively.  The  characteristic  indirect  gaps  are  clearly 
seen.  These  band  structures  compare  favorably  with  the  best  results  available,  and,  in 
contrast  with  those  obtained  in  the  usual  empirical  TB  approaches,  produce  good  con¬ 
duction  bands. 


III.  ALLOY  CALCULATION 


A.  VCA 

Because  we  have  the  same  Hq(1()  matrix  for  both  silicon  and  germanium,  it  is  only 
the  Hi  matrices  of  the  constituents  in  scaled  VCA  which  distinguishes  them.  In  this 
approximation,  the  diagonal  elements  of  the  alloy  Hamiltonian  H(k)  are  simply  the 
concentration-weighted  average  of  the  corresponding  elements  of  the  pure  silicon  and 
germanium  Hamiltonians,  whereas  the  off-diagonal  elements  of  R(ic)  are  obtained  by 
assuming  a  1/d^  dependence.  n(ic)  can  be  diagonalized  to  obtain  the  VCA  band  struc¬ 
ture  for  various  concentrations,  x.  The  VCA  band  structures  for  x  =  0.1  and  0.5  are 
shown  in  Fig.  1(c)  and  1(d)  respectively. 

B.  Diagonal  CPA 

An  earlier  work  on  CPA  band  structure  of  SiGe  alloys(^^)  is  based  on  a  local  but 
energy-dependent  pseudopotential  approximation.  While  the  value  of  the  scattering 
potential  parameter  was  1.49  eV,  close  to  our  value,  the  calculation  predicted  too-large 
linewidths  in  the  Eg  spectrum  and  essentially  no  effect  on  the  electron  mobility.  With 
the  availability  of  a  set  of  good  basis  functions  and  constituent  band  structures,  more 
realistic  band  structures  of  the  alloy  can  be  obtained. 

In  the  current  model,  we  have  a  TB  Hamiltonian,  which  contains  matrix  elements 
to  all  ranges.  The  simplest  alloy  model  is  to  assume  that  the  important  disorder  resides 
only  in  the  diagonal  matrix  elements,  f,  and  Cp.  In  our  model,  the  e,®'  and  e,®*  differ  by 


1.46  eV,  whereas  and  differ  by  0.21  eV.  For  the  present,  we  neglect  the  disorder 
in  the  off-diagonal  element.  Mathematically,  we  have 


(1) 


where  1  is  a  fee  lattice  vector  identifying  a  site,  and  Vy  is  the  8x8  diagonal  matrix  with 
elements  U,  =  c ,  -  T,,  Up  =  <p  -  7p  in  the  orthonormal  local  orbitals  |TjQf>;  j  denotes 
the  two  atoms  in  the  unit  cell  labeled  I,  a  represents  s  or  p  symmetry,  and  7,  and  7p  are 
the  concentration-weighted  average  values  of  s  and  p  silicon  and  germanium  term  value 
energies. 

The  one  particle  alloy  Creen’s  function  is  defined  as 

G«lloy(Z)  =  ■= — -  (5 


We  are  after  the  configuration  average  of  this  Green’s  function,  which,  in  effective 
medium  theory,  is  replaced  by  an  effective  Green’s  function  G, 


where  E  is  the  self  energy.  In  CPA,  we  can  now  write  E  =  (Ej),  with  Ej*  being  an 

T 

8x8  matrix  in  the  basis  |  lja>  having  the  form 


where 


E,  0  0  0 

0  Sp  0  0 

0  0  E.  0 

0  0  0  s. 


Here  E,  and  Ep  are  the  s  and  p  parts  of  the  seif  energy.  The  E,  and  Ep  are  determined 
from  the  conditions  that  the  average  atomic  t-matrix  with  respect  to  the  CPA  Green’s 
function  G  is  zero.  With  our  ansatz  for  E,  the  matrix  equation  <t>=0  reduces  to  two 
coupled  equations  <t,>  =  0  and  <tp>  —  0,  where  the  average  is  the  concentration- 
weighed  average  <Q>  =  xQ®'-t-yQ®*,  and  the  t  is  defined  as 

t/  =  (U/-EJ  [l-F„(U/-EJ]-‘ 


(a  =  s  or  p,  =  Si  or  Ge) 


In  the  above  expression,  is  the  diagonal  matrix  element  of  G  in  the  local  basis 
Fo(Z)  =  <Tjo  I  G(Z)  |Tjo>.  E,  and  Ep  are  coupled  because  F,  and  Fp  each  contain 
both  E,  and  Ep. 


An  iterative  average-t-matrix  (lATA)  procedure(*®)  is  employed  to  solve  the  CPA 
equation.  This  procedure  improves  upon  a  guessed  solution  through  the  following 


equations 


E.  =  E“  +  <t,‘‘>  [1  +  F“<t“>]-'  , 


Ep  =  Ep®  +  <tp‘'>  (1  +  Fp‘’<tp“>]-‘ 


(6) 


where  <ta>  and  F°  are  similar  to  those  in  Eq.  (5)  except  that  E^  now  replaces  E^. 
The  most  time-consuming  calculation  is  then  the  computation  of  the  local  Green’s  func¬ 
tions  Fg®  and  Fp ,  given  by  the  BZ  summation;  e.g. 


N  ?  lz-H{f)-E,®)„ 


where  the  inverse  of  an  8x8  matrix  is  involved  for  every  It.  This  can  be  simplified  by 
observing  that  E^  has  the  same  form  as  E  in  Eqs.  (4)  and  (5)  and  that  the  4x4  A  matrix 
can  be  written  as  A  =  Epl  (Ej-Ep)J,  where  I  is  the  identity  matrix  and 

-  ^ 


10  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 


Defining  the  matrix 


==  J  7) 


F,°{Z)and  Fp(Z)  can  now  be  calculated  from 


where 


W)  =  ^  E  «ii(^^.z)  , 


^  E  (822(1^, z)+g33(ir.z)+g«(ir,z)]  , 


g  =  g®  +  g'’(l-(7g°)-‘ 


8i°Ag.z)=E 


r  z-€„(k)-E; 


In  Eq.  (9),  €„(ir)  is  the  band  energy  in  VCA  and  {U,,  ^(Ic)}  satisfy  the  following  Eigen 
equation: 
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Because  the  a  matrix  has  only  two  nonzero  elements,  the  matrix  inversion  in  Eq.  (8)  is 
obtained  analytically. 

A  substantial  reduction  in  computer  time  is  made  possible  by  using  an  analytical 
continuation  method.(^^)  In  this  method,  £,  and  £p  are  calculated  as  a  function  com¬ 
plex  Z,  and  then,  using  the  analytical  properties  of  the  self  energy  and  Green’s  functions, 
they  are  interpolated  for  real  Z.  Because  the  functions  £„  £p,  and  G  are  smooth  for 
complex  Z,  the  CPA  iterations  and  BZ  integrations  can  be  carried  out  with  substantially 
less  computer  time. 

For  the  concentration  x  =  0.10  and  x  =  0.50,  the  L  and  X(A)  gap  respectively  are 
preferred.  The  L  to  X(A)  crossover  takes  place  near  x  0.15.  The  CPA  correction  to 
L  and  X  edges  at  x  =  0.10,  0.15,  and  0.50  should  be  good  enough  to  study  the  quantita¬ 
tive  variation  of  band  gap  in  Si^Gei.^  alloys.  Hence,  the  calculations  are  carried  out  for 
these  three  cases.  In  addition,  because  the  experimental  results  are  available  for  x  = 
0.109,  CPA  calculations  are  also  done  here  for  comparison.  As  expected,  E,  is  much 
larger  than  Ep  for  all  the  cases.  The  self  energy  as  a  function  of  energy  is  plotted  in 
Fig.  2  for  an  X  =  0.50  alloy. 


C.  Off-Diagonal  CPA 


As  mentioned  earlier,  silicon  and  germanium  differ  in  their  lattice  constant  by 
''•'4%.  In  order  to  include  the  effect  of  the  structural  disorder,  the  CPA  calculation  is 
repeated  next  with  off-diagonal  (OD)  disorder  included.  By  an  application  of  the  molec¬ 
ular  coherent  potential  approximation  (MCPA)(^^),  Hass  et  al.  included  OD  disorder  in 
the  CPA  calculation  of  A'  semiconductor  alloys.C^^)  Assuming  that  B  atoms 

occupy  the  sites  of  an  ordered  zinc-blende  virtual  lattice,  they  modeled  the  dominant 
structural  effect  as  the  difference  in  A'  -B  and  A'  '  -B  hopping  matrix  elements.  Hence, 
the  chemical  and  structural  disorder  effects  are  treated  as  random  variations  of  V|^, 
^2^,  where  the  symbols  have  their  usual  meaning.(35) 

The  extension  of  the  method  to  Si,Ge|_,(  alloys  is  not  straightforward,  mainly 
because  silicon  and  germanium  can  occupy  both  sublattices;  hence,  there  can  be  no 
ordered  virtual  lattice  in  this  case.  If  we  choose  the  tetrahedral  unit  cell  as  the  molecu¬ 
lar  unit  for  MCPA,  we  see  that  the  disorder  is  not  cell  diagonal.  However,  by  choosing 
an  appropriate  basis  set,  we  can  make  the  intercell  interaction  be  the  highest  order 
effect.  We  start  with  a  hybrid  basis  |Tli>j  obtained  from  the  sp®  hybrid  orbitals.(35) 
The  hybrids  1  through  4  (i  =  1-4)  are  obtained  from  orbitab  centered  on  a  sublattice  1- 
site,  and  the  states  5  through  8  (i  =  5-8)  are  those  from  the  orbitals  located  at  the  four 
nearest  neighbor  sites  on  sublattice  II.  The  Bloch  basis  states,  corresponding  to  Aj  T2 
symmetries,  |  k  >j  located  on  an  I-site  (i  =  1-4)  and  ll-site  (i  =  5-8)  are  obtained  from 
the  corresponding  hybrid  states  given  by  the  relation 


|T’>i=  E 

j=l 


(10) 


An  explicit  de6nitioD  o(  these  •ifbii  .  sn  be  found  id  Ref  (20)  In  this  new  basis,  the 
self  energy  E  at  the  given  «ite  takes  the  furm 


CTo' 


(1 


where 


<^o 


i:,  0  0  0 

0  0  0 

0  0  Ep  0 

0  0  0 


( 


E,'  =i(3E,  +  E.) 


(13) 


The  self  energies  can  be  obtained  again  from  the  lATA  iteration  procedure 


E  =  Eo  +  <<T>>  [1  +  F<<T>>  ]-‘  , 


B 


p  _  1  Q  ioQnj 

N  t  r  z-UiS) 


(R{k)  +  E)Q  =  E(Z,ir)Q 


<  <T>  >  =  X  <Ta>  +  y  <Tb> 


<Ta>  =  t*A,  +  t^A.B  +  *'*ab,  +  y^  t*B,  .  A  =  Si  ,  (16) 


with  a  similar  expression  for  <Tb>-  Physically,  for  a  given  A  atom  at  the  center,  the 
other  four  atoms  in  the  molecular  unit  cell  can  be  all  A  atoms,  3  A  atoms  and  1  B  atom, 
2  of  each  1  A  atom  and  3  B  atoms,  or  all  four  can  be  B  atoms.  <'r»>  or  <Tb> 
represents  the  conSguration-averaged  t-matrices,  and  <<T>>  is  the  concentration- 
weighted  average  of  the  configuration.  By  exploiting  the  symmetry,  as  seen  in  Eq.  (11), 
one  can  reduce  this  problem  to  solving  two  2x2  coupled  matrix  equations.  Eq.  (14)  can 
be  iterated  to  obtain  E„  Ep,  E^,  and  Ejp.  After  every  iteration,  we  get  a  new  set  of 
E,,  Ep,  Ejs,  E2p,  and  E|,'  :  The  new  set  has  not  been  tested  to  see  if  E^'  is  still  given 
by  Eq.  (13).  In  our  calculation,  we  did  not  iterate  to  obtain  a  new  Ej,'  ;  instead  we  fixed 
it  by  the  relation  given  in  Eq.  (13).  The  error  introduced  by  this  approximation  is 
expected  to  be  very  small.  As  in  the  case  of  diagonal  CPA,  the  computation  can  be  sub¬ 
stantially  reduced  by  the  method  of  analytical  continuation.(32) 


IV.  DISCUSSION 


A.  Eq  ar»d  Ej  Optical  Transition 

The  V'CA  values  of  Eg  (r2'  e  -  r2Sv)  and  Eq'  (ri5(  -  r25v)  and  their  measured 
values  are  plotted  as  a  function  of  x  in  Fig.  1(e).  Because  the  measurements(^®)  are 
made  at  room  temperature,  the  experimental  values  are  smaller  than  the  values 


.-a  <.(i_ 


calculated  from  the  zero  temperature  band  structure.  Inclusion  of  the  relativistic  effects, 
which  are  not  present  in  our  calculations,  is  expected  to  form  a  more  accurate  basis  for 
comparison  with  the  experiments.  As  seen  from  Fig.  1(e),  the  theoretical  and  the  experi¬ 
mental  values  both  have  a  linear  variation  with  x.  Similar  calculations  of  E|  (L|,  -  ,) 

also  have  a  linear  variation  on  x  and  are  in  qualitative  agreement  with  experiments.(^^) 

From  the  CPA  self  energies  E,  and  £p,  it  is  straightforward  to  calculate  the  correc¬ 
tion  to  the  VGA  bands.  The  calculated  complex  band  structure  is  plotted  for  x  =  0.50 
in  Fig.  (3).  The  CPA  corrections  are  shown  only  in  the  vicinity  of  the  band  gap.  The 
shaded  portion  represents  the  half-width  of  that  energy  state.  Because  s-scattering  is 
dominant  in  these  alloys,  we  see  that  the  major  disorder  lies  in  the  conduction  band. 

The  topmost  valence  band,  with  its  rich  p-content  is  least  affected.  The  CPA  band 
structure  is  used  to  calculate  the  Eq  and  Ej  peak  positions  for  x  =  0.10,  0.109,  0.15,  and 
0.50  concentrations.  The  calculations  and  the  data  from  Reference  (36)  show  a  small 
bowing  that  is  not  seen  on  the  scale  of  Fig.  1(e). 

The  self  energies  E„  Ep,  Ejp,  and  E^  are  calculated  in  MCPA  for  the  x  =  0.50 
alloy.  As  in  the  case  of  CPA,  the  self  energies  associated  with  s-symmetry  are  much 
larger  than  the  ones  associated  with  the  p-symmetry.  While  E25  is  found  to  be  very 
small,  Ejp  is  at  least  an  order  of  magnitude  smaller — almost  zero.  However,  the  ImE, 
obtained  by  CPA  and  MCPA  differ  considerably.  As  seen  from  Fig.  4,  the  difference 
increases  as  one  goes  away  from  the  band  edge.  Therefore,  the  lifetime  associated  with 
the  alloy  disorder  is  decreased  by  the  inclusion  of  OD  structural  disorder.  In  addition, 
the  OD  disorder  lowers  the  conduction  band,  introducing  an  extra  bowing.  The 
Eo  and  Ej  values  are  reduced  by  27  and  12  meV  respectively.  The  VCA,  CPA,  and 
MCPA  values  of  Eq  and  E,  are  listed  in  Table  III. 
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The  half-width  of  the  alloy  states  is  calculated  from  the  imaginary  part  of  the 
CPA  self  energies.  The  half-width  of  the  lowest-lying  conduction  band  of  Sio  jCcQs  alloy 
is  plotted  in  Fig.  5  as  a  function  of  in  the  [100|  direction.  The  calculated  half-width 
is  186  meV  for  the  r2'  state  and  decreases  to  zero  at  the  band  edge.  Because  of  the 
negligible  alloy  broadening  of  the  topmost  valence  band  state,  the  half-width  correspond¬ 
ing  to  the  Eg  transition,  A(Eo),  is  186  meV,  which  is  approximately  one-half  of  the  previ¬ 
ously  published  CPA  results.C^^)  The  CPA  value  of  the  half-width  corresponding  to  the 
El  transition,  A(Ei),  is  31  meV.  Because  of  the  increase  in  the  imaginary  part  of  the  self 
energies,  the  MCPA  values  of  the  half-widths  of  the  Eg  and  Ei  transitions  are  206  meV 
and  32  meV  respectively.  Because  the  complete  Eg  peak  is  not  shown  in  the  published 
electroreflectance  spectrum(^^),  it  is  difficult  to  estimate  the  corresponding  half-width. 
However,  one  can  conclude  from  the  spectrum  of  the  x  =  0.458  alloy  that  the  half-width 
of  the  El  transition  is  considerably  smaller  (  ^  50  meV)  than  that  of  the  Eg  transition. 
The  agreement  between  the  experimental  and  the  theoretical  values  can  be  regarded  as 
good  because  there  are  errors  in  estimating  the  width  from  the  published  spectra,  and  we 
have  neglected  the  extrinsic  broadening  due  to  the  apparatus  used  in  the  experiments. 

In  order  to  make  a  more  accurate  comparison  with  the  experiments,  the  CPA 
values  of  A(Eg)  and  A(Ei)  are  calculated  for  the  x  =  0.109  alloy.  The  calculated  half¬ 
widths  of  the  Eg  and  Ei  transitions  are  13  and  2  meV  respectively.  From  the  spectrum, 
we  estimate  the  corresponding  values  to  be  8  to  15  meV  and  3  to  6  meV.  We  see  that 
CPA  values  are  in  excellent  agreement  with  these  experiments.  Because  x  is  small,  the 
inclusion  of  off-diagonal  disorder  is  not  expected  to  change  the  calculated  values 
significantly. 
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B.  Energy  Gap 


The  fundamental  gaps  of  these  alloys  are  calculated  as  a  function  of  concentration. 
The  VGA  gap  is  an  increasing  function  of  x  with  a  slope  discontinuity  at  x  0.11.  The 
conduction  band  minimum  changes  from  L-point  to  X(A)-point  at  this  crossover.  In 
addition  to  the  band  gap,  the  effective  electron  masses  and  the  band  edge  Kq  are  also 
calculated.  When  the  X(A)  gap  is  preferred,  the  band  edge  moves  linearly  from 
k  at  {0.9,0,0),=(,  15  to  (0.8,0,0)x=i.  The  effective  masses  at  a  given  minimum  increase 
linearly  from  their  pure  germanium  values  to  the  corresponding  pure  silicon  values. 

Using  CPA  self  energies,  the  band  gap,  band  masses,  and  the  band  edge  are  also 
calculated.  The  position  of  the  band  minimum  did  not  change  by  virtue  of  the  inclusion 
of  off-diagonal  disorder.  While  the  effective  transverse  mass  remains  almost  the  same  as 
the  VGA  value,  the  longitudinal  mass  has  a  maximum  of  12%  enhancement.  Because 
the  real  part  of  CPA  self  energies  is  negative  in  the  forbidden  gap  region,  an  extra  bow¬ 
ing  is  introduced  to  the  VGA  energy  gap.  Because  of  this  bowing,  the  L-X(A)  crossover 
takes  place  near  x  ca  0.13.  The  VGA,  CPA,  and  experimental^®)  bowing  parameters 
are  0.06,  0.18  and  0.24  respectively.  The  calculated  energy  gap  is  plotted  as  a  function 
of  X  in  Fig.  6. 

Because  of  the  negligible  change  in  the  effective  masses,  the  corresponding  values  in 
the  pure  materials  are  used  in  the  calculation  of  the  alloy-scattering- limited  electron 
mobility.  The  CPA  X-gap  Eg^  and  L-gap  Eg^  are  fitted  to  a  polynomial  form.  The  gen¬ 
eralized  Brooks’  formula  that  is  applicable  to  the  alloys  with  an  indirect  gap  and  multi¬ 
ple  bands  is  uscd.(^®)  The  calculated  electron  drift  mobility  and  the  experimental  Hall 
mobility(*)  are  plotted  in  Fig.  7,  where  the  theory  explains  the  qualitative  behavior  of 
experimental  results.(*)  As  observed,^®)  even  a  few  percent  alloy  concentration  can 
reduce  the  drift  mobility  substantially.  It  can  be  seen  that  the  rate  of  decrease  near 


0  and  X  =  1  are  quite  different.  This  is  because  the  L-edge  has  more  s  content  than  the 
X-edge.  Because  the  s  scattering  is  dominant  in  these  alloys,  the  L  electrons  are  scat¬ 
tered  more  than  the  X  electrons.  Precisely  for  this  reason,  one  observes  a  dip  in  the 
mobility  near  the  L  to  X(A)  crossover.  For  x  <  0.13,  the  minimum  gap  is  the  L  gap. 
After  the  crossover,  the  minimum  gap  is  the  X(A)  gap,  and  the  reduced  alloy  scattering 
increases  the  average  mobility.  For  still  larger  x,  the  mobility  decreases  because  of  the 
increased  alloy  disorder.  All  these  features  are  clearly  seen  in  Fig.  7.  While  our  calcula¬ 
tions  include  the  inter-valley  scattering  mediated  by  alloy  disorder,  the  effect  of  other 
scattering  mechanisms  is  expected  to  increase  the  dip  near  the  crossover. 

The  calculated  alloy  scattering  rate  for  the  holes  is  several  orders  smaller  than  that 
for  the  electrons,  because  (1)  the  valence  band  edge  has  dominant  p  content,  (2)  the  p- 
scattering  parameter  (Afp  =  0.21)  is  only  1/7  of  Ae„  which  alone  decreases  the  scatter¬ 
ing  rate  for  holes  by  a  factor  of  50,  and  (3)  finally,  the  imaginary  part  of  the  self-energy 
is  proportional  to  the  density  of  states,  which  approaches  zero  at  the  band  edge.  Hence, 
the  hole  mobility  in  this  system  b  insensitive  to  alloy  disorder. 

In  MCPA,  the  conduction  band  is  pushed  down,  because  of  an  increase  in  the  ima¬ 
ginary  part  of  the  self  energy,  giving  rbe  to  an  additional  bowing  in  the  fundamental 
gap.  For  an  x  =  0.50  alloy,  the  gap  b  reduced  by  7  meV.  The  bowing  parameter, 
including  the  MCPA  correction,  b  0.21,  which  is  in  excellent  agreement  with 
experiment. 

It  is  interesting  to  compare  the  results  of  our  calculations  with  those  of  Hass  et 
al.(3‘^)  In  their  calculations  on  the  Ga,  Jnx.4s  alloy,  CPA  introduced  an  extra  bowing  in 
the  fundamental  gap.  However,  after  the  MCPA  corrections,  the  total  scattering  was 
diminbhed  and  the  results  were  similar  to  VCA  results.  These  results  were  explained  in 
terms  of  the  relative  strength  and  sign  of  the  atomic  term  values  and  We  extend 
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their  argument  to  Si^Gei.^  alloys.  The  hybrid  level  of  silicon  is  higher  than  that  of  ger¬ 
manium.  Because  of  its  shorter  bond  length,  the  V2  of  silicon  is  larger  than  that  of  ger¬ 
manium.  Thus,  in  this  case,  both  effects  combine  to  give  more  disorder  in  the  conduc¬ 
tion  and  valence  bands.  Therefore,  the  scattering  is  enhanced  in  these  alloys.  This 
explains  the  increase  in  the  imaginary  part  of  the  self  energy  due  to  inclusion  of  OD 
disorder  in  our  calculation. 


In  conclusion,  we  have  incorporated  both  chemical  and  structural  disorder  into  the 
calculation  of  the  CPA  band  structure  of  Sij^Gei.^  alloys.  The  calculation,  based  on  a 
realistic  band  structure  of  silicon  and  germanium,  suggests  that  the  band  gap  is  an 
increasing  function  of  x  with  a  slope  discontinuity  at  x  ^  0.13.  The  linewidths  of  the 
Efl  and  Ej  transitions  calculated  by  CPA  and  MCPA  are  in  good  agreement  with  experi¬ 
ments.  Addition  of  the  structural  disorder  to  the  diagonal  CPA  decreases  the  band  gap 
slightly  but  increases  the  s-part  of  the  self  energy  considerably  over  certain  energy 
ranges.  The  calculated  alloy-scattering-limited  electron-drift  mobility  is  in  qualitative 
agreement  with  the  observed  Hall  mobilities. 
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FIG.  1. 


FIG.  2. 


FIG.  3. 


FIG.  4. 


FIG.  5. 


FIG.  6. 


Calculated  VCA  band  structures  of  (a)  silicon,  (b)  germaniunn,  (c)  Sig.sGeo  ^ 
alloys,  and  (d)  Sio  joGeogo-  (e)  Calculated  VCA  values(solid)  and  the  experi¬ 
mental  values  (dashed)  of  the  Eq  and  Eg'  peak  positions  are  plotted  as  a 
function  of  alloy  concentration  x. 

The  variation  of  the  imaginary  part  of  the  seif  energy  E,  (CPA)  as  a  func¬ 
tion  of  energy  for  x  =  0.50. 

Calculated  CPA  complex  band  structure  of  the  Sio.6Geos  alloy.  Only  the 
bands  in  the  vicinity  of  the  energy  gap  are  shown.  The  shaded  portion 
represents  the  alloy  broadening. 

Imaginary  part  of  E,  (CPA)  and  E,  (MCPA)  as  a  function  of  energy  for  x  — 
0.50  alloy. 

The  variation  in  the  width  of  the  lowest  lying  conduction  band  as  a  function 
of  in  the  [100]  direction  for  the  x  =  0.50  alloy. 

The  variation  of  the  VCA  energy  gap  (dash-dotted)  and  the  CPA  energy 
gap  (solid)  as  a  function  of  x. 


FIG.  7. 


Calculated  drift  mobility  (solid)  and  the  experimental  (dashed)  Hall  mobility 
(Reference  1)  as  a  function  of  x. 
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